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SUMMARY 


A theoretical analysis is presented yielding sets of partial differential 
equations for determination of turbulent aerodynamic flowfields in the vicinity 
of an airfoil trailing edge. A four-phase interaction algorithm is derived to 
facilitate the analysis. Following input, the first computational phase is an 
elementary viscous-corrected two-dimensional potential flow solution yielding 
an estimate of the inviscid-flow induced pressure distribution. Phase C 
involves solution of the turbulent two-dimensional boundary layer equations up 
to the trailing edge, with transition to a two-dimensional parabolic Navier- 
Stokes equation system describing the near-wake merging of the upper and lower 
surface boundary layers. An iteration provides refinement of the potential 
flow induced pressure coupling to the viscous flow solutions. The final phase, 
if desired, is a complete two-dimensional Navier-Stokes analysis of the wake 
flow in the immediate vicinity of trailing edge. A finite element numerical 
algorithm is presented which is applicable to solution of all partial differen- 
tial equation sets of the inyiscid- viscous aerodynamic interaction algorithm. 
Numerical results of prediction in the wake of a NACA 63-012 airfoil are 
compared to experimental data. 



INTRODUCTION 


Turbulent boundary layer flows departing the trailing edge of an aero- 
dynamic surface are experimental ly verified to be strong sources of noise. 
Alteration of the aerodynamic surface can alter far^field intensity by local 
absorbtion as well as modifications to the turbulent flow structure prior to 
departing the trailing edge. A considerable volume of experimental data has 
been obtained for the case of a jet flow directed over a surface. Attention 
is now turning to aeroacoustic analysis of the near-field wake associated with 
airfoil flowfields induced solely by forward flight. 

This report presents a theoretical analysis yielding sets of partial 
differential equations governing turbulent aerodynamic flowfields in the 
vicinity of a trailing edge. It is a generalization of the parabolic 
procedure developed by Baker and Manhardt (ref. 1) for analysis of slot jet 
flows and downstream of the sharp trailing edge terminus of a planar flap. The 
solution procedure developed herein is a viscous-inviscid interaction algo- 
rithm and is applicable to analysis of subsonic two-dimensional airfoils at 
angle of attack. The primary requirement of this algorithm is accurate 
numerical prediction of the detailed distributions of mean and fluctuating 
velocity fields in the immediate trailing edge vicinity. 

To accomplish this goal, the four phase interaction algorithm presented 
in Table 1 has been developed. The input phase serves its standard purpose. 
Phase B is an inexpensive iterative sequence wherein the Laplacian (2D<|)) is 
repeatedly solved for the potential flow over an airfoil which is sequentially 
augmented in thickness distribution by computed estimates of the viscous 
boundary layer displacement thickness distribution 6*. This solution provides 
the initial estimate of the inviscid flow pressure distribution imposed upon 
the subsequent solutions. 

In Phase C, the consequentially more complete (and expensive) two- 
dimensional boundary layer (2DBL) and parabolic Navier-Stokes (2DPNS) equation 
systems are solved in the trailing edge vicinity. The 2DPNS solution yields 
detailed distributions of the mean and fluctuating velocity correlations, as 
well as the computed distribution of efflux velocity U 2 (xi) from the 2DPNS 
domain. The 2Dcfi solution from Phase B is repeated, using this "onset" velocity 
boundary condition specification, producing a refined inviscid flow pressure 
field. This sequence of 2DBL, 2DPNS and 2D^ solutions is repeated as necessary 
to yield a common solution for the coupling pressure distribution. 

The terminal Phase D solution is probably not required for sharp trailing 
edge airfoils at small angles of attack. It is required to determine the 
immediate vicinity trailing edge solution for a blunt-based airfoil. The 
requirement of Phase D is to refine the prediction of the 2DPNS solution within 
the immediate trailing edge vicinity, and is accomplished by the more expensive 
solution of the complete time-averaged, two-dimensional Navier-Stokes (2DNS) 
equations. An appropriate 2DNS solution domain is determined and boundary 
conditions applied thereon from the 2DPNS solution. The 2DPNS solution also 
provides the initial conditions for the iterative solution to 2DNS. 
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Table 1 



VISCOUS-INVISCID AERODYNAMIC INTERACTION ALGORITHM 


PHASE 

OPERATION 


A. 

Input 

1. Problem specified by Moo, C, t/C, a. Re 

E (Mach No., chord, thickness distribution, 
angle of attack, Reynolds No.) 

B. 

Vi scous-Corrected 
Potential 
Computati on 
(2D<),) 

1. Solve v2(j) = 0 on' basic airfoil and wake centerline 

2. Solve integral boundary layer equation for 6* 

3. Extrapolate 6* onto wake 

4. Repeat B. 1-3 until 6*(x) is stationary 

C. 

Turbulent 
Boundary Layer, 
Parabolic 
Solutions 
(2DBL/2DPNS) 

» 

1. Initialize Ui(x 2 ) from data or Cole's Law 

across 6(xi) at Xi/C =■ 0.9 

2. Establish polar coordinate system on 

0.9 < Xj/C £ 1.0 if required. 

3. Execute 2DBL on 0.9 < Xj/C < 1.0 using inviscid Cp. 

4. Initialize 2DPNS at xi/C = 1.0 

5. Execute 2DPNS on 1.0 < x^/C < 1.2 using inviscid Cp. 

6. Modify 2D(|) boundary conditions using 2DBL/2DPNS 

solution for U 2 "onset" velocity 

7. Compute 2D(j) to refine inviscid Cp distribution. 

8. Repeat C.2-7 until Cp stationary. 

D. 

‘ Complete 

i Navier-Stokes 
j Solution 

1 (2DNS) 

i 

1. Identify appropriate 2DNS solution domain spanning 

Xi/C = 1.0 + e 

2. Obtain initial -boundary conditions from 2DBL/2DPNS 

solutions 

i 3. Solve complete Navier-Stokes Equations 
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This report presents derivation of the appropriate differential equation 
systems including closure for turbulence. A finite element numerical algorithm 
is developed that is applicable to solution of each system of the identified 
differential equation systems. Numerical predictions with the algorithm, for 
Phases A through C, are presented and compared with experimental data for a 
NACA 63-012 airfoil at zero angle of attack. 


SYMBOLS 

a^ boundary condition coefficient 

c coefficient 

C chord length 

skin friction coefficient 
Cp pressure coefficient 

d ordinary differential 

0 diffusion coefficient 

e finite element index 

f function of known argument; coordinate curve 

g function; source term 

h metric coefficient; mesh parameter; integration step size 

H boundary layer shape factor 

k turbulence kinetic energy; polynomial degree 

K generalized diffusion coefficient 

s, differential operator; length scale 

L differential operator 

M Mach number; number of finite elements spanning 

n unit normal vector; nodes per element; dimensionality 
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N finite element interpolation function 

p pressure; generalized parameter; iteration index 

q generalized dependent variable 

Q generalized discretized dependent variable 

R domain of elliptic operator 

Re Reynolds number 

5 finite element assembly operator 

t time 

T temperature 

u-j velocity vector 

u, U reference velocity 

x-j Cartesian coordinate system 

y+ friction velocity Reynolds number 

a parameter 

Y ratio of specific heats 

3R boundary of solution domain r'^ 

6 Kronecker delta; boundary layer thickness; increment 

6* boundary layer displacement thickness 

A increment 

e turbulence dissipation function; convergence criteria 
e boundary layer momentum thickness 

transformed coordinate system 
5^. transformed coordinate system 

K Karman coefficient 

A multiplier 


V kinematic viscosity; general diffusion coefficient 

p density 

CT-. mean flow Stokes stress tensor 

* \i 

T Reynolds stress tensor; wall shear; differential element 

(j) velocity potential function 

X generalized initial -value operator 

ip conservation potential function 

finite element natural coordinate system 
0) turbulence damping factor 

global solution domain 

Superscripts : 

e effective value 

n dimension of R 

t turbulent 

T matrix transpose 

+ turbulence correlation function 

mass-weighted time-average 

time average 

" unit vector 

" mass-weighted fluctuating component; ordinary derivative 

Subscripts: 

“ global reference condition 

e finite element domain 

ijjjk,£ tensor indices 
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I freestream reference condition 

n normal 

0 initial condition 

t time derivative 

w wall reference condition 

Notation: 

column matrix 
square matrix 
union 

intersection 
belongs to 
absolute value 


{ } 

[] 

u 

n 


PROBLEM DESCRIPTIONS 


Overvi ew 


The several differential equation systems required for the viscous- 
inviscid aerodynamic interaction algorithm are each derived from the complete 
Navier-Stokes equations. The distinguishing feature of aerodynamic flows is 
their fundamental unidirectionality as induced by the forward flight of the 
aircraft. Of course, local regions of separated, i.e. omnidirectional, flows 
can be induced by adverse pressure gradients and/or abrupt changes in geometry. 
An example of the former occurs on the trailing edge upper surface of an air- 
foil at sufficiently large angle of attack. A blunt trailing edge illustrates 
the latter. All but the complete Navier-Stokes equation descriptions have 
imbedded in their derivation the assumption of unidirectionality. This 
assumption greatly simplifies the attendant numerical solutions as well as 
producing considerable cost savings in computer time. The interaction algo- 
rithm concept is basically an overt attempt to enjoy these benefits to the 
fullest extent. The appropriate differential equation systems are developed 
in this section. 


Time-Averaged Navier-Stokes Equations 

A state point in single species fluid mechanics is expressed within 
solution of a coupled nonlinear partial differential equation system describing 
conservation of mass, momentum and energy. Unique solutions can be obtained 
upon closure of this system by specification of constitutive behavior and 
boundary conditions. In Cartesian tensor notation the non-dimensional, conser- 
vation form of the Navier-Stokes equation system for isoenergetic flow is 


■•(o) 




( 1 ) 


L(pu. ) 


It 


3 

3X. 




pUjU. 


+ P6.j - Re' 


•1 

'^ij 


= 0 


( 2 ) 


The dependent variables in equations (l)-(2) have their usual interpretation 
in fluid mechanics, i.e. p is mass density, uj is the velocity vector, p is 
the static pressure, and a-jj is the Stokes stress tensor 


a • • = u 
IJ 


3X. 


3U_. 


3X. 


3 9Xj^ TJ 


(3) 
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where y is the dynamic viscosity. The equation of state for a perfect gas 
closes the definition, and the Reynolds number Re = p„‘U<„£/y^ where A is a 
characteristic scale length. 

Equations (l)-(3) are valid descriptions for both laminar and turbulent 
flows. For the latter, however, the solution becomes tractible in a practical 
sense only after time-averaging. While more definitive resolutions are the 
subject of current research, (c.f. (ref. 2), the conventional mass-weighted 
time-averaging is assumed adequate for the present requirements. Therefore, 
c.f. Cebeci and Smith (ref. 3), the Reynolds decomposition is defined as 


Uj(x.,t) = Uj(x. ) + Uj(x^. ,t) (4) 

The mass-weighted time-average velocity is defined as 

ijj = ^/p (5) 


and 


pu 


i 


1 im 



- pu^Odt = 0 


( 6 ) 


This operation yields the important relation 


pu.Uj = pu.Uj + pu^uj 


(7) 


Substituting equations (4)-(5) into (l)-(3), time averaging and collecting 
terms yields the time-averaged Navier-Stokes equations for isoenergetic flow. 


9(puJ 9 r. 1 


= 0 


( 8 ) 

(9) 


The solution of various simplified forms of equations (8)-(9) is required 
for the aerodynamic wake flow analysis. The following sections establish the 
appropriate sub-systems. 
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Viscous-Corrected Potential Flow 

The first requirement of the interaction solution algorithm is to 
compute a viscous-corrected two-dimensional potential flow distribution. 

Phase B in Table 1. The solution domain is a "sufficiently" large region 
surrounding the airfoil and the wake trajectory. The airfoil will become 
augmented in thickness distribution during the sequence of viscous corrections, 
see Figure 1. 



Ucc 



Figure 1- Illustration of Viscous-Corrected Potential 
Flow Solution Requirements 


The differential equations governing the two-dimensional potential and 
integral viscous boundary layer flows are trivial subclasses of equations 
(8)-(9). Define the gradient of the perturbation potential function as 

the difference between the reference and local velocity vectors, i.e. ^ 


U^. (Xj) E U„ 


9i{i 

ax. 


( 10 ) 


For all Mach numbers below transonic, and for all irrotational , isentropic, 
inviscid, steady and turbulence-free flows, equations (8)-(9) collapse to 
the elementary Laplacian 


L((l>) - = v2(j) = 0 

axt 


( 11 ) 
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The boundary condition for solution of equation (11) is established 
from the defining equation (10) by the inner product with the unit normal 
h^- as 


z(<l,) 


= 1±. 


9X. 


"i - 


"i"i 


ii-l 

Uoo U^.n^. 


( 12 ) 


If the contour with normal n^ corresponds to an inviscid streamline, for 
example, the boundary layer displacement thickness 6*, the last term in 
equation (12) vanishes identically. If not, the last term corresponds to 
the inviscid entrainment velocity component u^. In either instance, equations 
(11)-(12) represent a well -posed Neumann problem, and (j) = 0 is set at one 
location (only) on the 2Dcj) boundary. 

Segments of the potential flow solution domain coincide with the viscous 
displacement surface 6* and its extrapolation into the wake. The distribution 
of 6* is estimated during Phase B by solution of the integral form of the 
steady flow continuity and Ui momentum equations (8)-(9). The resultant 
differential equation is written on 9, the boundary layer momentum thickness, 
as (c.f . , ref. 3) 

In equation (13), u(xi) is the potential velocity at the boundary layer edge, 

H = 6*/e is the boundary layer shape factor, and Cf is the skin friction 
coefficient. Using the Ludweig-Tillman formula. 


e 0.246(10) (14) 

and the length scale e is used in evaluating the Reynolds number. 

The primary output of this computational sequence is the potential flow 
pressure distribution induced on the airfoil and wake flow solution domains. 




pu.u. 


W 


(15) 


Upon substitution of equation (10) and resolution into scalar components on 
the local normal-tangent coordinate system (n, s), with unit vectors (n, t), 
equation (15) becomes 


■ nV 


2|f (S.t) 



(16) 
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This completes definition of the equation systems required to be solved 
within Phase B of the interaction algorithm. 


Turbulence Closure Model 

The Phase C and D differential eq-ation systems require identification 
of a turbulence closure model, since time-averaging has introduced the dynamic 
Reynolds stress -pu^uj into the mean flow differential equation system (8)-(9). 
For the aerodynamic wake flows of interest, it is reasonable to assume 


pu;ui = -putuj (17) 

The exact differential equation for the kinematic Reynolds stress tensor 
u^u:j is (c.f. ref. 3, 4) 

' \J 


L(u7uj) - ^ (U5U<) + 


9U, 


3U. 


9X|^ ^-^k 9x^ 


2v 


9u:f.9u:: 

-f- 


au;: 

— - + 

9u:j’ 

J 

9X, 9X, 
k k 

p 


9X. 


9X| 


9uCu; 




= 0 


(18) 


where v is the kinematic viscosity y/p. An additional differential equation 
is required established, usually written on a turbulence dissipation rate 
funtion. Assuming dissipation occurs in the sub-scale and is isotropic, the 
corresponding turbulence dissipation rate e is defined by the contraction of 
the fourthterm in equation (18) as 


2 _ 9u:9ui 

3 ^ij^ = 9X^9X^ 


( 19 ) 
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The transport equation for dissipation function e is 


a. a ~ 

'-<=> ° n ^ 55^ <“k') * 


+ 2 


3^u:f ' 

2 

+ 9 


3X, 3X„ 

3X| 

'' ‘^k 3X^ 3X^ ■ p 3X^. 3X^. 

L k 

k 


= 0 


( 20 ) 


Equations (18) and (20) constitute seven partial differential equations, and 
are not closed since the third order correlations are undefined. Additional 
differential equations could be established, but they in turn would involve 
undefined higher order terms. Hence, modeling of third order correlations 
is usually invoked at a level of completeness dependent upon the dimension- 
ality and geometrical complexity of the physical system. 

Launder, Reece and Rodi (ref. 5) present a comprehensive analysis of a 
Reynold's stress closure hypothesis and numerical results for several allied 
classes of flows, see also Hanjalic and Launder (ref. 6). Gessner and Emery 
(ref. 7) report an order of magnitude analysis of the Launder closure model , 
for the parabolic flow approximation in Cartesian coordinates that yields an 
algebraic equation system for the Reynolds stress field in terms of dominant 
mean flow gradients and the turbulence kinetic energy field k. 


k e i u'.u'. (21) 

From tensor field theory, see also reference 8, the existence of a Reynold's 
stress constitutive equation is assured and of the form 


- u;ui - - «, + ••• 


( 22 ) 


Equation (22) displays the symmetry required of a stress tensor since 


Eik 


SXk 


(23) 


The a may be determined from some sub-dimensional f 1 ows . For example, in 
two-dimensional boundary layer flow wherein only - uju^ remains significant, 
the form of equation (22) must agree with the familiar two-equation closure 
model 


13 


(24) 


- Uiu£ 


= C 


k2 


V e 


3Ui 

8X2 


where C,j is an empirical constant (0.09). Comparing equations (22) and (24) 
yields «2 = C^k^/e. 

The tensor generalization of Gessner and Emery's analysis for three- 
dimensional parabolic flows can be established. The first three terms in the 
expansion of equation (22) are 


1 J 


ij 


J<i 

^ aOjI 

k3 

+ CoCi. V 

- 35 . ^ 

'3U. 

— - + 

3Uj 

£ 

3Xj 




3Xfc 


+ ... (25) 


This form is achieved using the model of Launder, et al (ref. 5) for the 
triple correlations in equations (18)-(20). The diagonal tensor a-- is 
defined as 




3k ^i^■j 


(26) 


The coefficients aj_ permit anisotropy in the turbulence normal stress field. 
For ai = Cl, a2 = 33 e C 3, the C-,- are defined by Launder, Reese and Rodi 
(ref. 5) as 

22(C^i - 1) - 6(4C^2 - 5) 

33(Cj;jl - 20 ^ 2 ) 

4{3C^2 - 1) 

C2 = 

iKCjrfi - 2C^2) 

22(C^i - 1) - 12(3C p - 1) 

33(Cj(i - 2C^2) 


44Cj(Si - 22C/iiC^2 - 128C^2 - 36C^2 + 10 

165(C^i - 2C^2)^ (27) 


Here C^i and Cjj2 are the "universal" empirical constants derived by Hanjalic 
and Launder (ref. 6). The suggested values are Cjjii = 2.8 and C^2 = 0.45. 


14 



The algebraic Reynold's stress model developed by Gessner and Emery 
(ref. 7) retained most first-order terms in the expansion. Equation (25), 
simplified in accordance with the parabolic assumption, but retaining all first 
and second-order terms yields 


1^2 


U 3 U 3 = 03 k - 2 C 4 ~ 


U 1 U 2 


U 1 U 3 = 


U 2 U 3 


- Cl 


E 

3x1 


3 U2 



e 

3x2 

isi 

3 U 3 

■ E 

3X3 


3 Ui 

e 

3X2 

iii 

3U2 


- Cu — 


E 3X3 

k 2 


3 U 2 3 U 3 

3 X 3 3 X 2 


^ k p k2 

Lo — Lil — ■ 

^ e ^ e 


''3U '' 


[_3X2j 


^3ti '' 


3X 


3 ; 


Co — Ch — 
^ e E 


3Ui 


9 x 2 


p k p k 2 
Co ”■ Cii 

E ^ E 


3Ui 


3x 


3j 


k^ 3Uj 3Ui 


C 2 g ^4 g. 3 X 3 


(28) 


The retained second-order terms are the C 4 -premultiplied expressions in uruj 

and ufuj. These terms provide the elliptic boundary value character required 
for solution of the transverse momentum equations as initial -boundary value 
problems. Under the parabolic order of magnitude analysis, Gessner et al. 
(ref. 7 ) obtained the parabolic form of the dissipation function as 


G = 


dUi 

^ + U{U| 


3 U 2 
3X^ 


(29) 


With equations (25) -(26) established, the aerodynamic interaction algo- 
rithm solution requirement for turbulence closure is considerably moderated, 
since only equation (18) contracted over i and j need be solved in concert with 
equation (20). Again employing the model of Launder, Reece and Rodi (ref. 5), 
the final form is 
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L(k) - 1^-+ ® 


at ax . 

J 


u.k - 

( k 

C, - ulu'. + v6. . 

] Sk 

J 

k e 1 J 1J 

J 3x. 


32 ,- 


( 30 ) 


I ^ \ _ a E , a 
'■(=) - at 


F “i'-j i, 


au. 2 

+ u'.u' f + C2 |— =0 

E 1 j k ax. e k 

<J 


(31) 


The various C|^ and C“ are empirical model constants. 


Parabolic Solution Algorithm 

Following completion of the Phase B viscous-corrected potential flow 
solution, the next requirement is solution of the two-dimensional parabolic 
Navier-Stokes (2DPNS) equations for the pre- and post- trailing edge flow- 
fields. The familiar boundary layer equations are a subset of the 2DPNS 
system. There are four assumptions required satisfied for 2DPNS to be a valid 
flowfield description, specifically: 

1. The flow is steady. 

2. A dominant mean flow direction is omnipresent. 

3. In this direction only, diffusion is negligible compared to all 
other processes including convection, and 

4. Overall elliptic boundary value character can be enforced by pressure 
coupling with a complete-dimensional inviscid (potential) flow 
solution. 

The steady viscous and turbulent aerodynamic flows under study meet all four 
requirements except in a region of separation as induced either by a sufficient 
adverse pressure gradient or abrupt change in geometry. Excluding these, the 
3DPNS partial differential equation system is readily established from the 
complete Navier-Stokes equations (8)-(9), and turbulence kinetic energy and 
dissipation equations (30)-(31) as 
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( 32 ) 



(33) 

(34) 

(35) 


(36) 


L(e) 



(pU^e) 



C3 uju,; 


k 

e 



1 3U. 

" f 3^ " = 0 (37) 

Equations (32)-(37) are generally valid on two- and three-dimensional 
space dependent solely upon summation index limits. The 2DPNS limited summa- 
tion index convention is 1 £ i £ 2 and = k = 2, while 3DPNS 1 £ i £ 3 and 
2 £ s,, k £ 3. The Xj axis is assumed aprallel to the predominant flow 
direction, and the various remain the model correlation coefficients. 
Furthermore, the dilitation term from the Stoke 's laminar stress tensor, 
equation (3) has been deleted from equations (33)-(35) as negligibly small. 


For the Ui equation, upon inserting and expanding terms in the Reynolds 
stress, the divergence term becomes 


3 

y 

3U 1 

3 

fp + c pkM 

3X„ 

Re 

3X„ ■ 


iRe ^ 0 3x„ 

a 

_ 

1 

£ 



(38) 


Defining an effective diffusion coefficient 



(39) 
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the diffusion expression for in 2DPNS is the familiar form 


3 


> 1 

3U2 


_ 3 
3X„ 

- eff 

3X 

M 

Re 


- PU 2 U' 

- 

1 


a 

L K 


( 40 ) 


The divergence terms in the remaining momentum equations are 


3 

p 1 


* “" 4 ) 


3 

_ eff 



Re 

[ax. 



^"2 , 

pv 



3 

3X. 


k k 

Pl<«|^ - PC2 7 - k — ^ 


(41) 


where addition of terms involving shear of the alternative transverse plane 
mean velocity component u^ and the dominant component u^ are now present. 

The divergence term common to both the k and e equations contains three 
terms. After inserting equation (2'5), the first term in the expansion is 


3 

C p — u"u" 

3 ' 

|c p 

3X. 

q^ E 1 £ 3X, 


3 e aXj, 


(42) 


where q represents either k or e. Note that equation (42) is fundamentally 
identical to the Reynold's stress contribution to equation (39) with C 4 
replaced by Cq. For 2Cq/3 ~ C 4 , equation (39) yields the effective viscosity 
as the representative major diffusion contribution throughout 2D PNS. When 
inserted into equation (36), the second term 
identically in the parabolic approximation, 
a non-vanishing dif fusion- type contribution 


in the expansion of ufuj vanishes 
However, the third term yields 
for both the k and e equations as 


3 

k k k 

C D - Co - Ch - k 

' 3 Ui' 

puil 


3 X^ 


[axj 

laxj 
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Hence, an anisotropic divergence-basedtransport of both k and e is permitted 
for a non-symmetric three-dimensional iii velocity field. 


The production terms for both k and e are basically identical, with the 
, fkl-i 

latter multiplied by — . The Reynold's stress-velocity shear contraction 

simplifies to the following form under the parabolic approximation. 




3 X,. 


pCit 


Jsi 

3 U 1 3 U 1 

e 

3 X„ 

_ £ £_ 


+ 2 pC 2 7 C 4 f k 


'aui' 

~ 3 Ui 3 Ui' 

3 Xi 

3 X^ 3 XJ, 


(43) 


The first term is identical to the familiar boundary layer form. The second 
term may be negligible due to the appearance of aui/3xi « 3ui/9x^. 


Pressure Resolution and Mass Conservation 


The procedure to enforce pressure coupling and computational conservation 
of mass is required. As described in Table 1, the first step following 
initiation of a Ui velocity profile in Phase C is integration of the 2DBL 
form of the 2DPNS equations. This procedure establishes an initial distribu- 
tion of U 2 that satisfies the continuity equation, as well as to initialize 
distributions of k, e and veff in the absence of data. In 2DBL, equations 
(33), (36) and (37) are solved for Ui, k and e, while the continuity equation 
(32) is solved directly for U 2 as 


L(p) 


3 (pU2 ) 


3 X: 


3 (pUi) 

+ r- = 0 


3 X 


1 


(44) 


Equation (44) is an initial-value problem, the direct solution of which 
specifically enforces continuity. The integration is initiated at the 
airfoil surface using the specified value of (suction) velocity Vw , i.e., 
U 2 (xi, 0) = VwUi)- The second, term in equation (44) is the axial (x^) 
derivative of u^, an accurate evaluation of which results following a few 
integration steps of equations (33), (36) and (37). 

The pressure distribution through the boundary layer in 2DBL is estab- 
lished from the familiar order of magnitude analysis of equation (34) as 


L(pU2) 


3 

3X2 


+ o 


p U 2 U 2 




(45) 


Equation (45) is a linear ordinary differential equation solved exactly in 
terms of the appropriate Reynolds stress component. 
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Within the 2DPNS equation system, pressure coupling to the freestream 
inviscid flow, and the distribution throughout the viscous flow domain, is 
contained within solution of a transverse plane Poisson equation. Forming 
the divergence of equations (34) plus (35) for steady flow, and recalling 
the limited summation convention, this equation is 


UP) = 4 




9 

9X, 






( 46 ) 


The aerodynamic interaction algorithm is based upon the observation that the 
solution of the linear, elliptic Poisson equation (46) consists of comple- 
mentary and particular solutions as 


p(x,) = P^(x,) + Pp(x,) 


(47) 


By definition, the complementary solution satisfies the homogeneous form of 
equation (46) 


UPc> = 



= 0 


(48) 


subject to the known boundary conditions for p(x.), which corresponds to the 
bounding inviscid flow pressure distribution. Tl'ierefore, on 2DPNS domain 
segments coincident with the inviscid flow, 

P(x^) = (49) 

where x^ indicates x^ constrained to the boundary of the 2DPNS solution domain. 
Elsewhere, the appropriate boundary condition is vanishing normal gradient. 

UPc)=4"*'° 

The particular solution Pp is any function satisfying equation (46) with 
homogeneous boundary conditions on 2DPNS boundary segments coincident with 
the freestream. Elsewhere, a gradient boundary condition can be established 
from the inner product of equations (34)-(35) with the surface normal. Since 
the convection term vanishes identically at a wall or symmetry plane, the 
constraint is 
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(51) 


3Pd ^ 
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3X, 




- pu. 


Since the laminar viscosity term is negligible everywhere away from a wall, 
equation for Pp can be simplified to the inviscid form. 


L(Pp) = 







= 0 


(52) 


Equations (47)-(52) define the pressure-coupled interaction algorithm, 
as well as a solution procedure for pressure distributions in the equation 
solutions. Since the current requirement is solely two-dimensional, the 
2DPNS algorithm becomes identical with 2DBL for flows upstream of the trailing 
edge terminus. Since a = 2 only, equation (48) becomes 



dx^ 

2 


0 


(53) 


Applying p^. = p^ at the freestream and equation (50) at the wall, the solution 
is P(, = constant = pi, i.e. the complementary pressure is the inviscid flow 
pressure imposed across the boundary layer. Similarly, since Pp vanishes at 
the freestream, the solution of equation should agree essentially with the 
boundary layer form, equation (45). 

This redundancy does not hold for two-dimensional flows downstream of the 
trailing edge terminus for non-zero angles of attack. On Xi/C > 1, the upper 
and lower surface inviscid pressures are equal only for the planar wake 
associated with zero angle of attack. Designating the levels as pi and pi, 

equation (48) becomes a two-point boundary value problem, the solution to 
which is a linear interpolation between the end point values. Hence, the 
distribution of p^ in 2DPNS is known a priori. Equation (52) is solved 
directly with Pn = 0 at each end point. The x^-derivatives in the non-homo- 
geniety cannot be deleted since local strong accelerations occur parallel to 
Xi immediately downstream of the trailing edge. 

The final requirement within 2DPNS is explicit computational enforcement 
of conservation of mass. An implicit enforcement occurs in the non-homo- 
geniety in the particular pressure Poisson equation by direct insertion of 
the continuity equation. The explicit enforcement is accomplished within the 
finite element solution algorithm by considering the continuity equation as 
a differential constraint in the PNS solution statements for transverse 
velocity. Looking to the variational calculus for guidance, a suitable measure 
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of mass conservation is a potential function ip, wherein the Laplacian of ip 
is equated to the error in exact satisfaction of the continuity equation (32), 
by the computed mean velocity mass flux, i,e. 


L(^) = 


3 / -.A/ \ 


( 54 ) 


The boundary conditions for p are arbitrary and may include the derivative 
constraint 




= 

3X^ £ 


= 0 


(55) 


as well as setting = 0 at select locations. In this regard, P exhibits 
much of the character of velocity potential function associated with the 
2Dp solution. However, in the context of a differential constraint, the 
gradient of P acts in the manner of a pressure gradient in the PNS solution 
statement. Importantly, the PNS solution is acceptable only upon equation 
(54) becoming computationally homogeneous, whereupon i/j = 0 throughout. Hence, 
interpretation of the gradient of ij; as a velocity or pressure contribution is 
correspondingly a computational zero. 

Finally, for algorithm stability in the computation of subsonic flows, 
the particular pressure distribution computed during the current iterate of 
the overall interaction algorithm is applied solely within the transverse 
momentum equation. The previous iterate distribution is added to the current 
distribution of complementary pressure to form the x^-pressure gradient for 
Ui. In the limit of convergence, the computed particular pressure distribu- 
tion takes on a computationally stationary value. 


Coordinate Systems 


Only for an airfoil with.cusped trailing edge can the locally Cartesian 
boundary layer coordinate system be carried directly into the wake. For a 
non-zero trailing edge included angle, it is necessary to transform to a 
polar coordinate system well upstream of the trailing edge terminus, to 
"rotate" the upper and lower boundary layers to a planar coincidence at the 
trailing edge. Figure 2 illustrates the geometry for a sharp edged airfoil, 
and Xi/C = 0.9 is selected as the point to initiate transformation to polar 
coordinates. The intersection of the two rays, perpendicular to the local 
tangent to the airfoil at x^/C = 0.9 and 1.0, defines the origin of the polar 
coordinate system. (Only the lower surface system is shown; the operations 
on the upper surface are identical.) The exterior potential flow solution. 
Phase B, or experimental data on ui determines the boundary layer thickness 
6 at Xi/C = 0.9. The 2DBL/2DPNS domain intersection with the lower free- 
stream is extrapolated on x^/C > 0.9 and into the downstream wake as 
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the curve noted, fj. The computational requirement is that the region defined 
by this curve, and the airfoil surface or wake trajectory, contain sCxi). The 
extrapolation can be refined, if required, during the iteration sequence within 
Phase C, but it is not required to be a coordinate surface of the polar coordi- 
nate system x^- = {e,r}. Specifically, this boundary is assumed expressed in 
the form 


fi(xi) = ri (1 + ai0) 


(56) 


where F'l is the radial coordinate of fi(xi/C = 0.9), and ai is an input para- 
meter to produce the desired scaling. The other surface of the 2DBL domain 
is the airfoil contour, assumed expressed as 


f2(xi) = f2(l + a20) 


(57) 


On the interval 0.9 £ Xi/C £ 1.0, the 2DBL/2DPNS equations are cast onto 
the Xj a { 0 ,r} coordinate system. The predominant velocity component is Uq. 
For example, equation (33) takes the form. 
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(58) 


Derivatives in the remaining members of each equation set are correspondingly 
expressed. A grid stretching coordinate transformation, that uniformly scales 
the viscous flow domain on the interval bounded by fi and f 2 , is 


Q' = Q 

r' = ~ ^i(e) 

[f 2 ( 0 ) - fi(e)]f^ (59) 


where f is a normalizing factor. Using the chain rule, the 0-derivative in 
equation (58) becomes 
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( 60 ) 
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where superscript prime denotes the e-ordinary derivative. Dividing- by r, and 
repeating the operation for the radial derivative, yields the final forms 


i 9_ 
r 30 


i 1_ 

r se" 



3 3, - u 3 

3r" 3r ~ ^i3r" 


(61) 


The h-j, 1 £ i ±3 in equation (61) are functions of the metric of the 
coordinate transformation, and are defined as 


ff/(f2 - fl) 1 




■{hjai 


hj(a2 - ai)/f 


(62) 


Therefore, the 2DBL/2DPNS equation sets are recast onto the fixed {e^,r''} 
coordinate system. For example, all members of both sets possess common 
convection and diffusion operators. Equation (58) projected onto {Q',r'} 
serves the illustration as 


L(u,) 


-r, 

r 30' 


- pu. 
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(63) 


The transformation-induced radial pressure gradient on Pc vanishes 
identically for 2DBL, but only in 2DPNS for zero angle of attack. Note also 
that the influence of the Xi-(ie., 0 ) gradient of the inviscid pressure p^. 
is diminished as r~i through the boundary layer. Therefore, even though 
Pq equals a constant in the 2DBL solution, its influence is distributed as 



r 30' 


dxi 



1 + r'/(hifi) 


(64) 
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where dpVdxj is the potential flow pressure gradient. 

The upper and lower surface 2DBL solutions are computed independently up 
to the trailing edge. Thereat, the 2DPNS solution domain is defined as the 

interior of the region bounded by fi and fi, see Figure 2. A return to 
rectangular Cartesian coordinates is permitted since the upper and lower 
surface boundary layers have acquired a common alignment. Alternatively, if 
the wake is curved, due to angle of attack for example, a single polar coordi- 
nate system would be useable. An alternative is to employ a grid-translation 
transformation that retains fine grid resolution in the region with extremal 
gradients. This transformation, which also permits grid stretching, is simply 
the Cartesian equivalent of equation (59), ie. 


X2 


Xl 

Xz - fl(Xi) 


[f2(Xi) - fi(Xi)]f^ 


(65) 


While the fj could be general curves, piecewise linear segments of the form 
axi + b should be adequate. Then, 


S _ 9 
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( 66 ) 


Elliptic Solution Algorithm 


Analysis of separated flows in the trailing edge vicinity usually require 
a complete two-dimensional Navier-Stokes (2DNS) solution to accurately resolve 
local velocity distributions. Analysis of the trailing edge region of a blunt- 
based airfoil at zero angle of attack is an example, and Figure 3 illustrates 
the geometry. The cylindrical coordinate system is the descriptor for the 
2DBL/2DPNS solutions that provide initialization of the 2DNS solutions. 

The switch-over from dual 2DBL solutions, to the single 2DPNS solution for 
the sharp trailing edge airfoil, occurs at x^/C = 1.0 - n, where n > 0 provides 
the computational flexibility for the transition. For the blunt-based airfoil, 
an inviscid "extension" is extrapolated to a sharp intersection, defined to 
occur at x^/C = 1.0 as shown in Figure 3. The 2DBL solution is pursued on the 
interval Xi/C ^0.9 up to the blunt base terminus. Each 2DPNS solution is 
initiated just prior to this location, and continued in the {0" r"} coordinate 
system up to Xi/C = 1.0. The upper and lower 2DPNS solutions are then coupled 
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in a rectangular coordinate description and marched into the wake. The 2DPNS 
upper and lower surface extension solutions are solved independently as is the 
case for the 2DBL solutions. The no-slip wall boundary condition on Ui is 
removed on the extension, allowing the individual 2DPNS high shear regions to 
accelerate prior to actual merging. 

The solution domain for a Phase D solution of the 2DNS equations lies in 
the immediate vicinity of the blunt-based trailing edge as illustrated in 
Figure 3. The two-dimensional Navier-Stokes equations for an isoenergetic 
flow are contained in equations (8)- (9) by the limited summation index conven- 
tion 1 £ i, j 2. The solution of the 2DNS equations for small Mach number 
and large Reynolds number is a difficult task. The 2DNS algorithm suggested 
is based upon the formulation developed for 2DPNS, extended to three-dimensional 
space, ie., 3DPNS. Specifically, equation (9) in expanded form is 

9u . 3u . 9u . 9p 

Lfu,) = 5 ^ + eu, ^ + 5D2 ^ 
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3Xi 



pU^Uf 



pu:u£ 


(67) 


The pressure field remains resolved according to equation (47). The solution 
for complementary pressure is given in equations (48)-(50) with 1 ^ £ 2. 

The particular pressure is the solution to equations (51)-(52) with all sub- 
scripts allowed to range (1, 2). The x^ derivatives contained within the 
j-summation in equation (52), for the 2DPNS formulation, are now the time 
derivatives from equation (67). The 2DNS velocity field uj is determined 
from application of the continuity equation as a differential continuity 
constraint. 

Equation (67), and the corresponding 2DNS k and e equations (30)-(31), 
restricted to two-dimensions, require Reynolds stress components for closure. 
The constitutive equation (25) is assumed a valid representation. In the 
blunt base region, all terms in equation (25) may be required retained, since 
the assumptions yielding the simplified forms in equation (28) are specifically 
parabol ic. 

The initial 2DNS solution domain is bounded by the inviscid extension 
employed for the 2DPNS solution and the base. Boundary conditions for 
{q} = {ui, U 2 , k, e, Pq, pp> are provided thereupon by the 2DPNS solution. 

Upon a satisfactory 2DNS approach to steady-state, the 2DNS solution domain 
should be enlarged, as illustrated in Figure 3. The flow conditions on the 
new boundary are provided by the previous converged 2DPNS solution, and the 
2DNS system is resolved on the interior of the expanded domain. The primary 
purpose in the final 2DNS solution is a more accurate prediction of the corner 
flow at the base. However, the discretization required to produce an adequate 
resolution is a significant limiting factor in terms of computer time and 
storage. 
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FINITE ELEMENT SOLUTION ALGORITHM 


The several sets of partial differential equations that are required 
solved for the developed aerodynamic interaction algorithm, have been identi- 
fied. For Phase B, equations (11)-(13) are appropriate. The 2DBL/2DPNS 
system for the parabolic Phase C operations are given by equations (32)-(37), 
(44), (48)-(52), and (54)-(55). The elliptic Phase D sequence requires 
solution of equations (67), (30) and (31) in addition to (48)-(52) and (54)- 
(55). Fortunately, in terms of developing a numerical solution procedure 
for the aerodynamic interaction algorithm, each member of the set belongs to 
the general class of second-order elliptic partial differential equations 
expressed 



+ 


fi(q 



P. X.) 


+ f2(q>x) 


= 0 


(68) 


In equation (68), q is any dependent variable, ie. {q} = {u^. , k, e, p^., Pp, ij>}. 
K is the generalized diffusion coefficient, f^ is a function of its argument 
that includes convection and any parameter p, and f 2 is an initial-value 
operator if present. For the 2DNS equations, x is time t and 1 i 2, 

For the 2DBL/2DPNS equation system, x corresponds to the Xj coordinate and 
s, = 2, 1 £ i ^ 2. For the solution variables pp, p^, and ({), f 2 vanishes 
identically in all instances. 

In all instances, the solution domain corresponds to the product 
of an n-dimensional Euclidean space spanned by the x^ coordinate system 

1 £ A _< n, with a one-dimensional space spanned by the x coordinate, ie. 


X X e x^ X Xq»x) (69) 


The boundary 3ft of the domain ft is its intersection with another enveloping 
solution domain, or the infinite freestream, and 


X„, x) (70) 

On 3ft the boundary conditions appropriate for solution of all variables are 
contained within the expression 


3ft = 3R X X e X^ X 
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(71) 
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Here, the a^*(x£,x) are specified coefficients appropriate to egch q. Except 
for p^, Pp and <j), an initial condition is required specified for each q as 


The finite element numerical solution algorithm for equations (68)-(72) 
assumes all q and p interpolated on sub-domains fig, which form the finite 
element discretization of fi, as 

qg(Xj^,x) = {N|<(x£)>^^Q(x)}e (73) 


Determination of the sum of the expansion coefficients {Q}p> which constitute 
the numerical solution, is accomplished using the Method of Weighted Residuals 
in the form. 
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{0} 


(74) 


where Sg is the assembly operator. The multiplier x is evaluated to enforce 
cancellation of generated surface integrals with the prescribed boundary 
condition normal gradient. For qg = u only, eq 0 and the last term expresses 
the explicit enforcement of the continuity equation (32), ie. equation (54) is 
a differential constraint on the transverse momentum equation solution. 

Inserting equation (68) into (74), and using the divergence theorem on the 
lead term, produces the equivalent matrix solution statement 
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The global rank of equation (.75) is identical with the total number of 
node points on RnUaR at which the dependent variable requires solution. 

For f£ non-vanishing, equation (75) is a system of first-order ordinary differ- 
ential equations. For f 2 vanishing identically, it is large-order algebraic 
and the matrix structure is symmetric, sparse and banded. In this instance, 
standard matrix solution procedures are directly applicable. 

An implicit numerical integration algorithm is suggested to solve equation 
(75) for the various initial-value problem descriptions. Completing the terms, 
the equivalent matrix statement for equation (75) is 






Se 

[C]^{Q}g + 

[Ul, [K]^_ 

tQ}e + {f>e 


{ 0 } 


(76) 


where the square brackets denote square matrices of element rank, the curly 
brackets denote column matrices, and the superscript prime indicates the 
ordinary derivative with respect to x- For the first three terms, which account 
for acceleration, convection and diffusion respectively, the explicit matrix 
forms are 


[Clg = 
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{N} {Nl'dx 


Xl 


[U], = {PU,>; 
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(77) 


All terms not explicitly involving the dependent variable {Q} are contained in 

'«e- 
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An accurate single-step implicit integration algorithm is the trapezoidal 

rule 

" I [«>j+l " «>$] (78) 


In equation (78), j is the initial -value coordinate index, and h is integration 
step-size Ax- Follov/ing the usual matrix manipulations, insertion of equation 
(76) into (78) yields a large order, non-linear algebraic equation system. The 
Newton matrix iterative algorithm for solution of this system is 



The dependent variable in equation (79) is the iteration vector, and 


( 80 ) 


v/here p is the iteration index. The right side of equation (79) is the 
homogeneous form of equation (78) evaluated with the pill iterate, i.e., 

(F}P,1 = Se[('le(«>J+l «>j) " l(‘9e>J+I * '9e>j)J (81) 


where 

^ + [K]J {Q)P + {f)^ ( 82 ) 


In equations (81)-{82), the matrix operations constitute inner products on 
matrices of element rank, with the assembly operator yielding the equivalent 
global expression. The vanishing of {F}to within definition of a computed 
zero, yields equation (79) homogeneous, hence convergence of the iteration for 

any evaluation of the Jacobian. The initial estimate (QJj+j each iteration 

is typically 
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By definition, the Jacobian is the derivative of equation (81) with 
respect to (Q}P. It can be computationally evaluated as 


[J] = 




[U], 


( 83 ) 


All operations again involve matrix inner products of an element rank. The 
rank of [J] equals the order of {6Q}; Dirichlet boundary constraints are applied 
within the evaluation of {F}. 


NUMERICAL RESULTS 


The developed numerical solution algorithm has been implemented into a 
computer program for the parabolic analysis for a sharp trailing edge airfoil. 
Phases A-C of the interaction algorithm. Table 1, have been assessed for 
accuracy and convergence for a NACA 63-012 airfoil with 0.61m chord for which 
experimental data are available for comparison (ref. 9). 


Viscous-Corrected Potential Computation 

The potential flow solution procedure is an extension of that reported in 
reference 10.. The basic macro-element, and resultant finite element discreti- 
zation of an appropriate region of R^ is illustrated in Figure 4, as well as 
delineation of gradient boundary conditions. The solution of the integral 
boundary layer equation (13) yields the distribution of momentum thickness 
e(x), hence displacement thickness 6*(x), which in the subsequent potential 
flow solution augments the airfoil thickness distribution t(x)/C. The integral 
boundary layer solution is terminated at the trailing edge. The wake distri- 
bution of 6* is a quadratic interpolation between the trailing edge value and 
zero thickness at 1.5C downstream. 

Table 2 summarizes the results of the typical first three iterations of a 
Phase B analysis, for the NACA 63-012 for a = 0° and C = 0.61m. The streamwise 
stations span 0.776 £ Xi/C £ 2.5, where x^/C = 1.0 is the trailing edge. The 
Phase B solutions for pressure coefficient converge to ± 0.0001 Cp after 
approximately seven iterations, and Table 3 summarizes typical results. At 
zero angle of attack, differences in Cp between natural laminar to turbulent 
transition, and tripping at 5% chord, are small for a = 0°. The 6* augmentation 
procedure is equally applicable to non-zero angle of attack; results obtained 
for a = 6° are illustrated, and convergence was achieved in seven iterations. 
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Table 2 


Viscous-Corrected Phase B Potential Flow 
Pressure Distributions, Iterations 1-3 
NACA 63-012 Airfoil, C = 0.61m, a = 0° 



Pressure Coefficient 

Cn(xi) 1 

5tati on 


Iteration Number 

xi/g 

1 

2 

3 

0.776 


-.0231 

-.0216 

0.841 


.0448 

.0466 

0.895 


.0985 

.1004 

0.940 

.1361 

.1391 

.1406 

0.975 

.1458 

.1639 

.1632 

1.000 

.1410 

.1783 

.1741 

1.050 

.0903 

.0846 

.0831 

1.180 

.0592 

.0560 

.0553 

1.390 

.0213 

.0214 

.0212 

1.680 

.0090 

.0088 

.0087 

2.050 

.0017 

.0017 

.0017 

2.500 

.0017 

.0016 

.0016 


Table 3 

Viscous-Corrected Phase B 
Potential Flow Pressure Distributions 

NACA 63-012 Airfoil, C = 0.61m 



PT'essure 

Coefficient Distribution, Cp(xi 

) 

Station 

Zero Degrees Angle of Attack 1 

Six Degrees Angle of Attack 

Xi/C 

Natural Transition 

Tripped at 5%C 

Upper Surface 

Lower Surface I 

0.776 

-.0217 




0.841 

.0463 




0.895 

.1001 




0.940 

.1402 




0.975 

.1632 

.1608 


.1789 

1.000 

.1746 

.1715 


.1575 

1.050 

.0833 

.0817 


.0868 

1.180 

.0557 

.0546 

.0683 

.0434 

1.390 

.0212 

.0209 

.0453 

-.0039 

1.680 

.0087 

.0085 

.0395 

-.0278 

2.050 

.0017 

.0017 

.0270 

-.0240 

2.500 

.0016 

,0016 

-.0049 

.0084 
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Low Turbulence Reynolds Number Model 

The primary requirement of the attached flow turbulent boundary layer 
analysis near the trailing edge is to establish initial conditions for the 
parabolic Navier-Stokes solutions in the wake. In turn, a suitable procedure 
is required to develop initial condition profiles for the boundary layer 
solution distributions. Cole's Law (ref. 3) is available to establish a Ui 
mean velocity initial profile, and the continuity equation (32) yields the 
corresponding U 2 distribution. In certain instances, experimental data may be 
available to initialize the turbulence kinetic energy initial profiles. If 
corresponding Reynolds shear stress data are available then equation (29) 
could be employed to initialize isotropic dissipation function. 

Several very practical problems emerge in establishing suitable initial- 
condition distributions for the 2DBL equations using the TKE closure procedure. 
Aerodynamic boundary layers are typically thin, hence experimental measurements 
of either mean or fluctuating velocity correlations are extremely difficult to 
obtain near the surface. For example, for the comparison data set summarized 
in the Appendix, hot wire probe dimensions limited the nearest approach to 
approximately 0.01m, at which location Ui/Ug = 0.7. Hence, most of the detail 
of the Ui and k distributions eludes the experimentalist, and alternative 
methods are required developed. Furthermore, and of paramount impact for the 
current study, the numerical solution domain is required to reach to the wall, 
since flow evolution in the wake centerline region is specifically of primary 
interest. 


A viable alternative approach involves use of mixing length theory to 
predict the near wall distributions of turbulent viscosity vt. As is well 
known, and recalling equation (39), 


vt = C4 


e 


( o )£)2 


3 U]^ 
3X2 


(84) 


where (co£) is the Van Driest-damped Prandtl mixing length z, and s. is a piece- 
wise linear function of X 2 . The functional form ofco, equation (84) is known 
and depends strictly on various specified parameters (ref. 3). Importantly, 
equation (84) depends solely on geometry and Uj , which is available from Cole's 
law in the absence of detailed data. For the current study, the TKE closure 
model was extended to the wall by assuming a low turbulence Reynold's number 
modification based upon Van Driest damping concepts. The ratio (k/e) was 
assumed unaffected, and k and C 4 were modified in equation (84) as k => w k 

and C4 ==> C4, where 0 < a £ 2 is a parameter to be optimized. Hence, 

the distribution of k in the near wall region, 0 < y+ < 80 could be modified 
with the resultant dis tribu tion of C 4 affecting primarily the k and e source 
terms, as dominated by uiua, see equations (36)-(37). 

A sequence of numerical experiments is continuing on evaluation of the 
concept. One particularly demanding test case corresponds to the Bradshaw 
relaxing flow boundary layer, IDENT 2400 of reference 11. Figure 5 summarizes 
the accuracy of the basic MLT solution in terms of familiar boundary layer 
integral parameters. The standard value of the MLT parameter X = 0.09 produced 
correct overall solution trends, but detailed discrepancies from data. Increas- 
ing to X = 0.11 produced the quite acceptable agreement indic ated^ Using 
A = 0.11 to compute the corresponding initial distribution of U 1 U 2 , the flow 
was recomputed using the TKE closure model, see Figure 6 . 
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Fig. 5. Prediction of Bradshaw Turbulent Boundary Layer Test Case 
Mixing Length Turbulence Closure Model. 











As observed, for a = 1, use of the "standard" TKE model constants Cg = 1.3, 

= 1.44, and = 1.92 produce essentially correct trends but detailed 

inaccuracy. Setting = 1.92 considerably improved the agreement as indicated. 

Figure 7 illustrates corresponding computed distributions of mean velocity and 
Reynolds shear stress at the initial and final stations, which are separated by 
a distance of 1.3m. Viewing the distribution of Ui , the near wall region under- 
goes significant acceleration and the final station solution is in excellent 
agree ment with data. Over the same interval, the distribution of shear stress 
ufu 2 flattens considerably. These distributions are smooth as computed from the 
Uj, U 2 , k and e solutions , usi ng the derived constitutive relation equation (28). 
The deviation of the sum \x\\x\H about k was bounded by two percent throughout 
the solution field. 


NACA 63-012 Interaction Analysis 

The Phase A-C interaction analysis has been exercised for a NACA 63-012 
airfoil for M„ = 0.09, a = 0° and C = 0.61m. A thoroughly comprehensive data 
set has been obtained for this configuration on 0.90 < xi/C £ 1.10, see 
reference 9. The Appendix lists these data at the first and last stations, 
and at two stations immediately bounding the trailing edge terminus. Figure 
8a) -8d) are composite summaries of the measured mean and fluctuating velocities 
on the interval 0.9 < Xi/C < 0.9979. All profiles are by and large parallel; 
the flattening of the HTuT distributions at the last two stations give indica- 
tion of anticipation of the imminent trailing edge terminus. 

The Phase C in teracti on a nalysis was initiated using the experimental data 
for Ui, u'l s (uiu^p, and UTuI on the given interval at Xi/C = 0.90. Following 
considerable experimentation. Cole's law, MLT concepts, and the turbulence con- 
stitutive equation (28) were successfully employed to fill in the initial pro- 
files between the wall and the first experimental data points. The 2DBL solution 
was then exercised on 0.9 £ Xi/C < 1.000, and the resultant free-stream distri- 
bution of U 2 inserted into the potential flow analysis. Phase B, to modify the 
inviscid flow Cp. Four sequences through this interaction, using 2DPNS to 
supplant 2DBL, were sufficient to render the freestream distribution of U 2 
stationary to within + 2%. 

Figure 9 illustrates the evolution of the final 2DPNS solution on 
0.9 £ x/C £ 1.0, on a plot scale directly comparable to the experimental data 
distributions. Figure 8. In Figure 9, however, each curve is shifted to the 
left to retain definition of near_wall g radie nts. Recall the data employed to 
initialize the solution were Ui, ui, and UiUa, with use of MLT concepts to 
complete the variable distributions. Figure 9a) illustrates the evolution of 
Ui, for the noted Xi/C stations. The experimental data at Xi/C = 0.9 and 0.9979 
are plotted for comparison on the first and final computed solutions. Agreement 
at the first station is essentially by definition, and the final station com- 
parison shows quite good agreement for the selection of a = 1 in the low turbu- 
lence Reynolds model and the dissipation equation coefficients = 1.67 and 

= 1.92. 
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(a) Mean Velocity uj 

Fig. 8, Experimental Data for NACA 63-012, a = 0° , Ref. 9. 
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(d) Reynolds Stress - ufuf 


Fig. 8, Concluded. 
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(a) Mean Velocity Uj 

Fig. 9 Computed 2DBL/2DPNS Trailing Edge Solution. 
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DOWNSTREAM VELOCITY 


In Figure 9b), the plot of uf, the agreement with data at the first 
station is not exact. However, the 2DPNS solution has actually marched two 
percent chord. The computed solution recovers some of the initial disparity 
and the final station comparison is essential parallel curves with computation 
nominally 5% lower than data. The elevation in uf adjacent to the wall per- 
sists throughout the solution. Figure 9c) shBws the comparison for u| at the 
selected Rations. No agreement at the first station is required since the 
data for u£ are not used for initialization. However, by the final station, 
the comparison between data and numerical prediction is acceptable^ The near- 
wall peak in computed Un is somewhat subdued compa red t o that for uf. A simi- 
lar and modest near-walT peak is also computed in ufU 2 . see Figure 9d. A 
modest disparity between data and prediction is illustrated at the initial 
profile and increases in magnitude and extent in proceeding to the final station 
where the maximum difference is about 20%. The magnitude of the near-wall peak 
remains nominally constant throughout the solution. Figure 9e is the plot of 
dissipation length scale 


^d 



(85) 


where ^ <5) is defined equal to the value at X 2 = 6. 

In the overall view, the agreement between prediction and experiment is 
favorable and the final station solution fields for Uj , U 2 » k and e should be 
acceptable as initial data for the 2DPNS wake solution. The near wall model 
and the dissipation equation correlation constants can exert a profound 
influence on the numerical solutions. Figure 10 summarizes the impact of 
variation of model constants on the distributions of fluctuating parameters 
at Xi/C = 1.0. The dashed curves correspond to the results of Figure 9 for 
which C| - = aCj, = 0.25. The solid curve was obtained for aC^, = 0.35, 

while tne dash-dot curve corresponds to aC^. = 0.20. The level of aC^, reflects 
the difference between annihilation and production of dissipation, see equation 
(37). The larger levels tend to raise sharply all stress levels near the wall 
and to flatten the profiles away from the wall (near freestream). 

Since the data case corresponds to a = 0°, these computed trailing edge 
lower boundary layer distributions can be reflected about the mean chord to 
initialize the 2DPNS wake solution on a discretization containing 80 elements. 
For solution comparison. Figures lla)-lld) are composite plots of the experimen- 
tal data distributions on the interval 1.0029 Xj/C 1.10. The dominant 

^tion in the field near the trailing edge appears manifested primarily in 

U 2 which exhibits a noticeable peak which decays rapidly. The profiles of uf 
and U 1 U 2 also exhibit peaks with significantly smaller amplitudes. 
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Fig. 10. Influence Of Dissipation Equation Constants On 

Computed Reynolds Stress Distributions, Xi/C = 1.0. 
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(a) Mean Velocity Ui 

Fig. 11, Experimental Data for NACA 63-012, a = 0°, Ref. 9- 
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Using freestream pressure distribution of the final 2DBL/2DPNS iteration, 
2DPNS wake solutions were generated on 1.0 < x^/C < 1.10. Figure 12 summari- 
zes the computed solution which was obtained by elTminating the low turbulence 
Reynolds influence at the wall and setting the turbulence dissipation correla- 
tion coefficients to their standard values, = 1.3, 1.44 and 1.92. In Figure 
12, the experimental data profiles have been over-plotted at the nearest compu- 
tational profile. Viewing the distributions of mean velocity Ui, good agree- 
ment occurs at Xi/C = 1.0037, but the 2DPNS solution progressively under 
predicts acceleration of the momentum deficit on the centerline. At Xi/C = 1.10, 
the maximum discrepancy is 8%, There also is evidence of modest under predic- 
tion of the lateral spread of the wake. 

This is more graphic in the solutions for the fluctuating velocities uT 
and TjT, Figures 12b)- 12c). For the former, recall that the trailing edge 
2DPNS solution under-predicted data by about 5%, see Figure 9b). The experi- 
mental data indicates a nominal 10% increase in the level of uT, uniformly over 
the profile, immediately downstream of the trailing edge (compare Run 55 curve 
in Figure 8b) to Run 25 curve on Figure 12b)). There is no mechanism in the 
numerical solution to elevate the initial data profile in this manner; hence, 
the nominal 20% difference persists th roug hout Figure 12b). Referring to the 
constitutive relation, equation (28), UiUi is dominately Cj k, hence the 
levels of turbulence kinetic energy are correspondingly low. This is also 
evident in the freestream comparison of UT , hence may account for the noted 
under-prediction of wake spreading. 

The computed distributions of Uf and "uT about the wake centerline do exhibit 
a significant difference that is confirmed by the data. Both fluctuating 
velocities peak immediately downstream; thereafter, TIT’ exhibits the character- 
istic double hump while "uT acfiieves a nominally flat central plateau. For the 
computations, the sole difference b|tween the two velocity distributions, see 
equation (28), is in the term 2C^ ~ . For two-dimensional incompressible 
flows, the continui ty e quation (32) indicates that 9ui/3xi = - 9u2/9x2. Hence, 

the second term in uTuT can be replaced by the negative of the second term in 

U 2 U 2 . The resultant differences in uT and U 2 are in general agreement with the 
data which tends to confirm a viability for the Reynolds stress constitutive 
equation. 

Figure 12d) compares the 2DPNS predictions for Reynolds shear stress TTfuJ 
to data. The low levels away from the center region can be attributed to overall 
under-prediction of the level of k. The development of the center region peak, 
and its decay and lateral spread, are in qualitative agreement with the data. 
Figure 12e) shows the computed distributions of turbulence kinetic energy, the 

profiles of which are essentially comparable to uf , 
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Fig. 12. Computed 2DPNS Wake Flow For NACA 63-012, a = 0° 











CONCLUSIONS 


A theoretical analysis has yielded an interaction algorithm for prediction 
of mean and fluctuating velocity distributions in the vicinity of the sharp 
trailing edge of a subsonic airfoil. Governing systems of partial differential 
equations have been established and a numerical solution algorithm identified. 

The results of numerical prediction for a NACA 63-012 airfoil have been 
compared with detailed experimental data. While some detailed quantitative 
differences exist, the results of the numerical solutions have exhibited quite 
notable qualitative agreement in delineation of differences in scalar compo- 
nents of the Reynolds stress tensor. Within the theoretical framework of the 
derived. constitutive equation, exhibited differences in the wake fluctuating 
velocity distributions have been correlated directly with the action of the 
continuity equation. Interestingly, the parabolic order of magnitude analysis 
indicates this term is higher order, yet its impact on correlation with data is 
most important. 

A considerable effort is required to generate initial data for the numeri- 
cal predictions. The computational solutions are quite sensitive to boundary 
data, and in particular, a key required element was establishment of a low 
turbulence Reynolds number model, to account for the influence of the airfoil 
surface on the differential equation solutions for turbulence kinetic energy 
and isotropic dissipation function. The solution field of primary interest 
appears essentially dominated by the action of this model, minor modifications 
to which can cause truly significant solution differences. It is quite apparent 
that considerable theoretical attention should be applied to this specific 
topic, to render the developed procedure more widely applicable. Certainly 
there is little desirability in developing the Phase D analysis until this has 
been accomplished. 
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APPENDIX 


Select Experimental Velocity Profiles For 
NACA 63-012 Airfoil, a = 0°, U» = 30 m/s, C = 0.61m 
Untripped Boundary Layer 


Run Numbers Chord Station, xi/c 


34,35 

0.9000 

54,55 

0.9979 

24,25 

1.0029 

22,23 

1.0052 

12,13 

1.1000 
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RUN MO.- 34 Y-STATION- 


1 PROBE- 2 VriAX- 0.99312E+O1 


POSITION (IN U UT 


UT 




6.77000E+«e 

e.au64E4-«6 

«.85328E-»ed 

0.89492E4^ea 

e.93£S£E+00 

e.978a«E+0e 

e.iei98E+ei 

0.ie615E+01 

•.11031E4^dl 

«.11448E^01 

e.ll864E-f91 

0.1S£8eE-f0I 

6.12697E+01 

e.l3113Em 

e.l3S3eE«91 

d. l39-<6E4ei 

e. l43SSE401 
6.14779E+01 
e.l5195E4ei 
9.i56iaE-»ei 
e.l£928E+01 
0.16444E401 
e.lS861E+dl 
0.17277E+ei 
e.l7694E4ei 
0.18110E401 
e.i9i5iE4ei 
e.2ei9SE4«i 
e.21233E401 
e.22274E+«l 
e.23315E401 
e.2435£E401 
0.2S397E401 
O.26438E401 
0.274798401 
0.285208401 
O.30602E401 
•.32684E401 
e.34?e6E401 
O.3S848E401 
O.38930E401 
0.41012E4O1 
0.43094E+O1 
0.4517SE401 
0.472SSE+01 
0.49340E401 
0.51422E401 
0.53504E401 
0.5553eE401 
0.576S8E4O1 
0.59750E4O1 
0. £39148401 
0.68078E401 
0.72242E401 
e,76406E401 
0.80570E401 
0.84734E401 
0.S8898E401 
0.93082E4O1 
0.97226E4O1 
0.10139E4O2 
0.105558402 
0.10972E402 
0.11S05E402 
0.12637E402 
0.13470E402 
0.14303E402 
0.15136E402 
0.159£9E402 


0.679168400 

0.682418400 

0.68565E400 

O.6S889E4O0 

0.692128400 

0.695358400 

0.69858E400 

0.701808400 

0.705028400 

0.708238400 

0.711438400 

0.714628400 

0.717808400 

0.720978400 

0.724138400 

0.727288400 

0.730418400 

0.733538400 

0.736648400 

0.739738400 

0.742818400 

0.745868400 

0.748908400 

0.751938400 

0.754938400 

0.757918400 

0.76529E400 

0.772538400 

0.779638400 

0.786598400 

0.793408400 

0.800068400 

0.806568400 

0.812918400 

0.819098400 

0.825128400 

0.836718400 

0.847708400 

0.858098400 

0.867948400 

0.877268400 

0.886118400 

0.894508400 

0.902488400 

0.910058400 

O.S1725840O 

0.924088400 

0.930548400 

0.936658400 

0.942408400 

0.947808400 

O.957S0840O 

0.965768400 

0.972638400 

0.978228400 

0.982668400 

0.986138400 

0.988798400 

0.990808400 

0.992328400 

0.993488400 

0.994408400 

0.995188400 

0.996628400 

0.997988400 

0.999088400 

0.999828400 

0.100078401 

0.100228401 


0.722448-01 
0.722728-01 
0.723008-01 
0.723278-01 
0.723528-01 
0.723768-01 
0.723988-01 
0.72418E-01 
0.724368-01 
0.724508-01 
0.724628-01 
0.724698-01 
0.724738-01 
0.724738-01 
0.724698-01 
0.724608-01 
0. 724468-01 
0.724278-01 
0.724048-01 
0.723758-01 
0.723408-01 
0.723018-01 
0.722568-01 
0.722058-01 
0.721488-01 
0.720868-01 
0.719078-01 
0.716928-01 
0.714438-01 
0.711608-01 
0.708468-01 
0.705018-01 
0.701288-01 
0.697298-01 
0. 693068-01 
0.688618-01 
0.679208-01 
0.669258-01 
0.658908-01 
0.648158-01 
0.636798-01 
0.624438-01 
0.610878-01 
0.595768-01 
0.579208-01 
0.561478-01 
0.542918-01 
0.523688-01 
0.503668-01 
0.482468-01 
0.459708-01 
0.409518-01 
0.357268-01 
0.309618-01 
0.268588-01 
0.231798-01 
0.197578-01 
0.166718-01 
0.141208-01 
0.122928-01 
0.111168-01 
0.101888-01 
0.920518-02 
0.771908-02 
0.707768-02 
0.653818-02 
0.610448-02 
0.564078-02 


0.628158-01 
0.626308-01 
0.624458-01 
0.622618-01 
0. 620798-01 
0.618998-01 
0.617218-01 
0.615478-01 
0.613778-01 
0. 612108-01 
0.610488-61 
0.608898-01 
0.607358-01 
0.605858-01 
0.604398-01 
0.602988-01 
0.601608-01 
0.600268-01 
0.598948-01 
0.597668-01 
0.596408-01 
0.595168-01 
0.593938-01 
0.592708-01 
0.591488-01 
0.590268-01 
0.587158-01 
0.583928-01 
0.580508-01 
0.576868-01 
0.572988-01 
0.568878-01 
0.564568-01 
0.560078-01 
0.555458-01 
0.550758-01 
0.541188-01 
6. 531498-01 
0.521588-01 
0.511268-01 
0.500348-01 
0.488798-01 
0. 476688-01 
0.464218-01 
0.451508-01 
0.438568-01 
0. 425278-81 
0.411378-81 
0.396648-01 
0. 380968-01 
0.36441E-01 
0.329888-01 
0.295928-01 
0.264508-01 
0.235548-01 
0.207988-01 
0.181168-01 
0.155528-01 
0.132598-01 
0.114068-01 
0.100158-01 
0.389578-02 
0.781728-02 
0.585568-02 
0.467728-02 
0.379258-02 
0.322508-02 
0.287948-02 


0.113918-02 
0.1 13828-02 
0.113728-02 
0.113628-02 
0.113518-02 
0.113398-02 
0.113258-02 
0.113118-02 
0.112958-02 
0.112778-02 
0.112588-02 
0.112378-02 
0.112138-02 
0.111878-02 
0.111598-02 
0.111298-02 
0.110968-02 
0.110818-02 
0.110248-02 
0.109348-02 
0.109428-02 
0. 108988-02 
0.108518-02 
0.108038-02 
0.107528-02 
0.107008-02 
0.105628-02 
0.104168-02 
0.102638-02 
0.101058-02 
0.994348-03 
0.977858-03 
0.951058-03 
0.943888-03 
0.926258-03 
0.908038-03 
0.869378-03 
0.827198-03 
0.781598-03 
0.733468-03 
0.684168-03 
0.634948-03 
0.586518-03 
0.538988-03 
0.492098-03 
0.445678-03 
0.399328-03 
0.355468-03 
0. 313078-03 
0.P7343E-03 
0.236908-03 
0.173278-03 
0.121528-03 
0. SI 4788-04 
0. 526748-04 
0.329588-04 
0.190238-04 
0.816318-05 
-0.316638-06 
-0.582768-05 
-0.839228-05 
-0.891178-05 
-0.833678-05 
-0.692718-05 
-0.611478-05 
-0.575608-05 
-0.533438-05 
-0.509918-05 


RUN NO.- 35 


TKE 


0.916508-02 
0.914588-02 
0. 912668-02 
0.910758-02 
0.908868-02 
0.906988-02 
0.905118-02 
0.903258-02 
0.901418-02 
0.899578-02 
0.897758-02 
0.895938-02 
0.894118-02 
0.892298-02 
0.890468-02 
0.888628-02 
0.886768-02 
0.884388-02 
0.882968-02 
0.881018-02 
0.879018-02 
0.876958-02 
0.874848-02 
0.872658-02 
0.870398-02 
0.868058-02 
0. 361818-02 
0.854938-02 
0. 847398-02 
0.839158-02 
0.830228-02 
0. 820668-02 
0.810528-02 
0.799898-02 
0.788868-02 
0.777518-02 
0.754188-02 
0.730378-02 
0.706208-02 
0.681488-02 
0.655858-02 
0.628898-02 
0.600398-02 
0.570418-02 
0.539328-02 
0.507598-02 
0.475618-02 
0.443478-02 
0.410998-02 
0.377898-02 
0.344128-02 
0.276528-02 
0. 215208-02 
0.165828-02 
0.1 27628-02 
0.959838-03 
0.718548-03 
0.519798-03 
0.375188-03 
0.281188-03 
0.223868-03 
0.182938-03 
0.145848-03 
0.938718-04 
0. 719688-04 
0. 571308-04 
0.476648-04 
0.401098-04 
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RUM HO.- S4 Y-STATIOH- 11 PROBE- 2 YflAX- O.92511E+01 
POSITION MN U UT UT U*U 


-0. 109102461 
e.ll3S6E401 
0.11743E401 
O.1S15SE401 
O.ia576E401 
e.ia993E401 
e.l3408E401 
O.1382SE401 
0.14S41E401 
d.l4658E401 
0.1S07-1E401 
0.154902401 
0.15907E401 
O.16323E401 
0.16740E401 
O.17156E401 
O.17572E401 
O.17989E401 
O.18405E401 
O.18822E401 
0.19238E401 
O.19654E401 
0.20071E401 
0.S0487E401 
0.209042401 
0.21320E401 
0.22153E401 
0.22986E401 
O.23818E401 
0,246512401 
0.254842401 
0.263172401 
0.271502401 
0.27982E4O1 
0.288152401 
0.296482401 
0.304812401 
0.313142401 
O.32979E401 
0.34645E401 
0.363102401 
0.379762401 
0.39642E401 
0.413072401 
O.42973E401 
0.446332401 
0.463042401 
0.479702401 
0.496352401 
0.513012401 
0.529662401 
0.54632E401 
0.562982401 
0.579632401 
0.596292401 
0.612942401 
0.644172401 
0.675402401 
0.706632401 
0.737862401 
0.789092401 
0.800322401 
0.83155E401 
0.862782401 
0.894012401 
0.925242401 
O.95647E401 
0.9877OE401 


0.716052400 

0.718352400 

0.720642400 

0.722942400 

0.725232400 

0.727522400 

0.729812400 

0.73210E400 

0.73438E400 

0.736652400 

0.738932400 

O.74119E40O 

0.74345E400 

O.7457OE400’ 

O.74795E40O 

0.75018E4O0 

0.752412400 

0.75463E4OO 

O.75684E400 

0.75904E4OO 

0.76124E400 

0.76342E400 

0.765592400 

0.76774E400- 

0.769892400 

0.772032400 

0.776262400 

0.780442400 

0.784572400 

0.788642400 

0.792652400^ 

0.796612400 

0.80051 2400 

O.80435E400 

0.808132400 

0.81185E400 

0.815S1E400 

0.81912E400 

0.826162400 

0.832982400 

0.839592400 

0.846012400 

0.852252400 

0.858322400 

0.864262400 

0.870052400 

0,875732400 

0.881292400 

0.886742400 

0.892092400 

0.897342400 

a.90849E4e« 

0.907522400 

0.91845E400 

0.917242400 

0.921912400 

0.930272400 

0.938052400 

0.945222400 

0.951732400 

0.957582400 

0.962782400 

0.967392400 

0.971472400 

0.975082400 

0.978292400 

0.981172400 

0.983772400 


0.682402-01 
0.681712-01 
0.681022-01 
0.680342-01 
.0.679662-01 
0.678992-01' 
0.678342-01 
0.677702-01 
0.677082-01 
0.676482-01 
0.675902-01 
0.675342-01 
,0.674802-01 
0.674292-01 
0.673812-01 
0.673342-01 
0.672912-01 
0.672492-01 
O. 672102-01 
0.671732-01 
0.671382-01 
0.671052-01 
0.670732-01 
-0.670422-01 
0.670112-01 
0.669812-01 
0.669202-01 
0.668542-01 
0.667782-01 
0.666892-01 
-0.665812-01 
0.664512-01 
0.662952-01 
0.661112-01 
0.658952-01 
0.656482-01 
0.653702-01 
0.650622-01 
-0.643642-01 
0.635792-01 
0.627312-01 
0.618432-01 
0.609232-01 
-0.599712-01 
0.589802-01 
0.679422-01 
0.568602-01 
0. 557512-01 
0.546382-01 
_ 0.535492-01 
0.625012-01 
0.51494E-01 
0.505062-01 
0.494982-01 
0.484262-01 
0.472502-01 
_C. 447202-01 
0.418532-01 
0.388382-01 
0.357612-01 
0.325882-01 
. 0.293612-01 
0.262852-01 
0.235642-01 
0.211782-01 
0.189262-01 
0.166942-01 
0.146262-01 


0.483892-01 

0.483562-01 

0.483222-01 

0.482882-01 

0.482552-01 

0.482232-01 

0.481912-01 

0.481602-01 

0.481292-01 

0.481002-01 

0.480712-01 

0.480432-01 

0.480162-01 

0.479902-01 

0.479652-01 

0.479412-01 

0.479172-01 

0.478942-01 

0.478712-01 

0.478482-01 

0.478262-01 

0,478032-01 

0.477792-01 

0.477552-01 

0.477302-01 

0.477042-01 

0.476452-01 

0.475782-01 

0.474932-01 

0.474052-01 

0.472962-01 

0.471692-01 

0.470252-01 

0.468622-01 

0.466812-01 

0.464832-01 

0.462712-01 

0.460462-01 

0. 455672-01 

0.450672-01 

0.445602-01 

0.440492-01 

0.435272-01 

0.429772-01 

0.423802-01 

0.417222-01 

0.410022-01 

0.402322-01 

0.394362-01 

0.386402-01 

0.378672-01 

O.3?126E-0l 

•.364132-01 

0.357072-01 

0.349862-01 

0.342292-01 

0.326822-01 

0.310092-01 

0.292882-01 

0.275332-01 

0.257092-01 

0.238502-01 

0.220872-01 

0.205142-01 

0.190692-01 

0.175972-01 

0.160542-01 

0.145812-01 


0.114882-02 
0.114902-02 
0.114932-02 
0.114952-02 
0.114992-02 
-0.115042-02 
0.115102-02 
0.115172-02 
0.115272-02 
0.115382-02 
O.US51E-02 
0.115652-02 
0.115822-02 
0.116012-02 
0.116222-02 
0.116452-02 
0.116692-02 
e. 116952-02 
0.117222-02 
0.117502-02 
0.117792-02 
0.118082-02 
0.118372-02 
0.118662-02 
0.118942-02 
0.119222-02 
0.119722-02 
0.120132-02 
0.120422-02 
0.120582-02 
0.120562-02 
0.120372-02 
0.120002-02 
0.119432-02 
0.118682-02 
0.117752-02 
0.116662-02 
0.115432-02 
0.112642-02 
0.109542-02 
0.106292-02 
0.102982-02 
0.996172-03 
0.SS183E-03 
0.926182-03 
0.838762-03 
0.849492-03 
0.808792-03 
0.767522-03 
0.726722-03 
0.687292-03 
0.849762-03 
0.614112-03 
0.579852-03 
0.546232-03 
0.512532-03 
0.447832-03 
0.382292-03 
0.319582-03 
0.262262-03 
0.210312-03 
0.163352-03 
0.122692-03 
0.903012-04 
0.662822-04 
0.483392-04 
0.339332-04 
0.222532-04 


RUN NO.- 55 


TKE 


0.699822-02 
0.698552-02 
0.697292-02 
0.696032-02 
0.694792-02 
0.693582-02 
0.692382-02 
0.691212-02 
0.690082-02 
0.688982-02 
0.637922-02 
0.686902-02 
0.685922-02 
0.684982-02 
0.684082-02 
0.683222-02 
0.682412-02 
0.681632-02 
0.680882-02 
0.680172-02 
0.679482-02 
0.678812-02 
0.678162-02 
0.677522-02 
0. 676872-02 
0.676212-02 
0.674832-02 
0.673302-02 
0.671542-02 
0.669462-02 
0,666992-02 
0.664072-02 
0.660642-02 
0.656662-02 
0.652132-02 
0.647042-02 
0.641432-02 
0.635332-02 
0.621912-02 
0.607332-02 
0.592082-02 
0.576492-02 
0.560632-02 
0. 544362-02 
0.527472-02 
0.509802-02 
0.491432-02 
0.472682-02 
0.454052-02 
0.436052-02 
0.419022-02 
O.40300E-O2 
0.387672-02 
0.372512-02 
0.356912-02 
0.340422-02 
0.306802-02 
0.271332-02 
0.236622-02 
0.203692-02 
0.172292-02 
0.143092-02 
0.117882-02 
0.976102-03 
0.812142-03 
0.667862-03 
0.536432-03 
0. 426542-03 


62 



RUN NO.- S4 


POSITION m 


».ioaasE-o3 

«.4t7<EE-01 

0.833S2E-01 

0.125«2E-^«» 

• .166SS£-f09 

• .20830000 
0.2-4904E400 
0.29158E400 
0.33322E+O0 
0.3748eE4^80 
0.41G50E4O0 
•.45S14E+O0 
0.49978E^00 
0.54142E400 
0.S8306E400 
•.e2470E«00 
0.72880E«00 
0.83290E+00 
0.937O0E«0O 
0.10411E«01 
0.11452E401 
0.12493£-»01 
0.13534E491 
0.14575Em 
0.15616E+01 
0. 17S98E-<-01 
0.19780E«01 
0.21862E401 
O.23944E401 
•.26026E«0t 
O.28108E+01 
0.30190EtOl 
O.32272E401 
0.34354E4O1 
O.36436E401 
0.38518E401 
0.40600E4O1 
t.42G82£4^01 
0.4e84££401 
0.Sl010Em 
O.55174E-^01 

• .59338E-»01 
0.63502Et01 
O.676G8E401 
O.71830E401 
O.7S994E401 
O.80158E401 
0.84322E401 
0.8848GE401 
0.9fiai4E401 
0.I0514E402 

9.11347C402 
0.1218OE402 
0.13013E402 
0.13845E4O2 
0.I4&78E402 
0.15SHE+02 
0.1G344E+02 
0.17177E+02 
0.180O9E402 
0.18842E402 
0. 19£75E402 


V-STATION* 


U 


0.4837GE400 

0.47674E4OO 

0.48969E40O 

0.50256E400 

0.51533E400 

O.5279GE40O 

0.S4O40E40O 

0.55262E40O 

0.58458E400 

0.57623E400 

O.S8755E400 

0.59848E40O 

O.60899E400 

0.S1905E400 

0.62863E400 

0.63769E400 

0.GS80OE400 

0.67487E400 

0.68847E400 

0.69924E40e 

0.70783E400 

0.71498E400 

0.72141E400 

0.727G5E4O0 

0.73403E+00 

0.74730E400 

0.76018E400 

O.77149E409 

0.78ie3E400 

0.79178E+00 

0.80238Etee 

0.81272E+00 

0.82208E400 

0.83O82E409 

0.84002E400 

0.85025C4Oe 

0.86092E4e0 

0.87094E409 

0.88727E400 

0.89982E490 

O.90861E40O 

0.91697E400 

0.9272OE40O 

0.93582E4OO 

O.94291E400 

O.95035E409 

0.95872E400 

0.9S654E40O 

0.97069E400 

0.98146E400 

O.98509E400 

0.g905SE400 

0.9945SE4OO 

0.99335E400 

0.99031E400 

0.98835E40O 

0.10035E401 

0.10055E401 

0.1O013E401 

0.10067E401 

0.10097E401 

O.1O086E401 


12 PROBE- 


UT 


0.71479E-01 

0.72655E-01 

O.73801E-O1 

0.7488SE-01 

0.75884E-01 

0.76771E-01 

0.77526E-O1 

0.78136E-01 

•.78591E>01 

0.78887E-01 

0.79028E-01 

0.79022E-01 

0.78884E-O1 

0.78G32E~01 

0.78290E-O1 

O.77882E-01 

0.76741E-01 

0.75709E-01 

0.75008E-01 

0.74669E-01 

0.74572E-01 

O.74536E-01 

e.7443eE-et 

0.74221E-01 

0.73964E-01 

0.73570E-01 

0.73299E-01 

0.72758E-01 

0.71994E-01 

0.71324E-01 

0.70782E-01 

0.70195E-01 

0.69437E--0t 

0.£85£1E~01 

0.67715E-01 

0.66798E-01 

0.65521E-O1 

0.63992E-01 

0.61744E-01 

0.S9844E-01 

0.57138E-01 

O.53288E-01 

0.49238E-01 

O.46212E-01 

O.41491E-01 

0.37117E-01 

0.33102E-01 

0.28678E-01 

•.24506E-01 

O.17489E-01 

0.13052E-01 

0.1022SE-01 

0.80144E-02 

0.64932E-02 

0.5S891E-02 

O.S4490E-02 

0.50132E-02 

0.46971E-O2 

0.46658E-02 

0.45649E-02 

0.45572E-02 

0.41435E-02 


2 YMAX- 


VT 


0.63700E-01 

0.64144E-01 

0.64527E-01 

0.64789E-01 

0.64879E-01 

0.64757E>01 

0.S4394E-01 

0.G3778E-01 

0.62907E-01 

0.61796E-O1 

0.60473E-01 

0.S8975E-01 

O.57347E-01 

0.55S40E-01 

O.53904E-01 

0.52190E-01 

•.4828OE>01 

0.45270E-01 

0.43317E-01 

0.42275E-01 

0.41836E-O1 

0.41685E-01 

e.41£0a£-01 

0.41489E-01 

0.41334E-01 

0.40993E-O1 

O.40689E-01 

0.4O363E-01 

0.40010E-O1 

0.39e24E-01 

0.39254E'01 

0.38908E-O1 

0.38496E-81 

O.38034E-01 

0.37634E-O1 

0.372e££-01 

0.36532E-01 

0.35B71E-01 

O.34283E-01 

0.32975E-01 

0.3t292E-01 

0.2g495E-ei 

0.27615E-O1 

0.25941E-O1 

0.23737E-01 

0.218O8E-O1 

0.2OO02E-O1 

0.1SO18E-O1 

O.16199E-01 

0.1304SE-01 

•.1044e£-01 

0.85708E-02 

0.6617OE-02 

•.5e8g8£>02 

0.41133E-02 

0.35733E-02 

0.31564E-02 

O.27429E-02 

0.24711E-82 

0.234S5E-02 

0.21SG7E-02 

0.30592E-02 


0.96S88E401 


U8U 


0.53782E-03 

0.24029E-03 

-0.48O63E-O4 

-0.318GOC-03 

-0.5637&E-03 

-0.77747E-O3 

-0.95547E-O3 

-0.1O955E>e2 

-0.11973E-02 

-0.12G25E-02 

-0.12944E-O2 

-0.12977E-02 

-0.12779E-O2 

-0.12410E-02 

-0.11932E-02 

-0.U401E-O2 

-01O14GE-02 

-0.9350GE-03 

-0.91148E-03 

-0.92614E-03 

-0.953£8E-03 

-0.97674E-03 

-0.99103E-03 

-0.10010E-02 

-0.10124E-02 

-0.10423E-02 

-0.I0G16E-92 

-0.10598E-02 

-0.10480E-02 

-0.10389E-02 

-0.10384E-02 

-0.19355E-02 

-0.10158E-02 

-0.987£4E-03 

-0.96618E-O3 

-0.94GG5E-03 

-0.91759E-03 

-e.88163E-03 

-e.81823E-03 

-0.74770E-03 

-0.G7226E-03 

-0.57848E-03 

>0.48£39E-03 

-0.41767E-03 

-4.32G49E-03 

HJ.25379E-03 

-•.19597E-03 

'O.14O01E-03 

-9.10013E-03 

-e.45223E>04 

-0.19182E-04 

-0.70776E-05 

-•.G7578E-0G 

0.2H51E-05 

0.21738E-05 

0.23232E-05 

0.18381E-05 

0.18178E-05 

0.17335E-05 

0.18892E-05 

6.18805E-05 

0.15894E-0G 


RUN NO.- 85 


THE 


0.91G69E-02 

0.93933E-02 

0.9G102E-02 

0.98055E-02 

0.9967GE-02 

0.10087E-01 

0.10J57E-01 

0.10173E-81 

•.iO134E-01 

0.10042E-01 

O.99084E-02 

0.97225E-02 

0.95114E-02 

0.92788E-02 

0.90350E-e2 

O.87S93E-02 

0.82201E-02 

0.77812E-02 

0.7G025E-02 

0.73G27E-02 

0.731 12E-02 

0.72932E-02 

0.727aSE-02 

0.72302E-02 

0.71793E-02 

0.70930E-02 

0.70276E-02 

0.69230E-02 

0.G7839E-02 

0.GG572E-02 

0.65511E-02 

0.G4411E-02 

0.G3034E-02 

0.61471E-02 

e.G001CE~02 

0.584G2E-02 

0.5G276E-02 

0.53674E-02 

0.4987GE-02 

0.4GG8GE-02 

0.42439E-02 

0.3709GE-02 

0.31870E-02 

0.2808GE-02 

•.22850E-02 

0.18532E-02 

0.14958E-02 

0.11471E-02 

0.8629EE-03 

0.47G15E-03 

0.27947E-03 

0.17801E-03 

0.10803E-03 

0.G80G8E-04 

0.51609E-04 

0.424G0E-'04 

0.35096E-04 

0.29S86E-04 

0.2787GE-04 

0. 26354 E-04 

0.25549E-04 

0.21409E-04 


RUN NO.- 2S V-STATION- . 13 PROBE- 2 VHAX- 0.9586SE+O1 «tf HO.- 23 
POSITION UN U UT VT U<U TKE 


•.336S2E-01 

O.7S262E-01 

•.11699E+0O 

0.tS854E-i-O0 

0.20018C400 

0.24182C+00 

O.28346E^00 

0.32510em 

0.36674E400 

0.40833E400 

0.45002E400 

0.49i££E-f00 

0.53330E400 

0.574S4E4^00 

t.61S58E+00 

0.65822E-^00 

0.6998SE400 

0.74150E-»00 

0.78314E400 

0.82478E-f00 

0.8£642E«00 

0.90806E-»00 

0.94970E400 

0.10538E^01 

0.11579E401 

0.12620E-^01 

0.13£61E-f01 

0.14702E«01 

0.15743C401 

0.16784E401 

f.l7825E«01 

0.18866E40I 

0.20948E401 

0.23030E^01 

0.2S112C401 

0.27194E401 

0.2927£E-»01 

0.313S8E401 

0.33440E-»0i 

0.35S22£-»01 

0.37604E+01 

•.3968SE«01 

0.417e8E40i 

t.43850E-»01 

0.45932E-^1 

0.S0096E-»01 

0.S42£0E-»01 

0.68424E401 

0.S2588E401 

> .SC7S2e<»0i 
.709tfiE4ti 

0.75080C401 

0.79244E401 

0.83408E401 

0.87S72E401 

0.9i736E401 

0.1000SE402 

0.10839E402 

0.11S72E402 

0.12505E402 

0.13338E402 

O.14170E402 

0.15093E402 

0.15838E408 

0.16669E402 

0.17502E402 

0.18334E402 

0.19167E402 

0.20000E+02 


0.57268E400 

O.57S93E-t-0O 

0.57923E400 

0.582C2E400 

0.586HE400 

0.58981E400 

0.593S8E400 

0.59776E400 

0.60207E400 

0.606filE400 

0.61139E400 

0.61638E400 

0.62157E400 

0.62S94E400 

0.63245E400 

0.63807C400 

0.64375E400 

0.84948E+00 

0.S5514E4O0 

0.8S075E40O 

0.S6827E40O 

0.S7163E400 

0.87682E400 

0.68885E400 

0.69934E400 

O.78827E400 

0.71581E400 

£.7223eE400 

O.72812E4O0 

0.73362E400 

9.73908E400 

0.744e4E400 

0.75621E400 

0.76797E4O0 

0.77939C400 

0.79O3iE>O0 

0.80086E400 

0.8n27E400 

0.82189E400 

0.83208E400 

0.84208E400 

0.8S115E-^00 

0.8S899E400 

0.86S91E4M 

0.87269E40O 

0.88744E400 

0.O0049E«80 

0.909S8E400 

0.91790C400 

t.9840ae4M 

0.93184E400 

0.942&8E400 

0.95321E40# 

0.96226E400 

O.97094E400 

0.97806E400 

0.988S5E400 

0.99442E40O 

8.99540E400 

0.99358E400 

•.O9640E400 

0.99878E4O0 

O.99871E400 

0.99992E400 

0.10014E401 

0.10015E401 

O.10019E4O1 

0.10O73E401 

0.10026E-^01 


0.64440E-01 

0.64215C-01 

0.S4029E-01 

0.C3920E-01 

•.C3920C-01 

0.fi40S3E-#l 

0.84337C*il 

t.fi4777E-01 

0.fi5372£-01 

0.66109C-01 

0.C6968C-01 

0.67923E-01 

0.68942E-01 

0.69991E-01 

f.71037E-Ol 

0.72047E-01 

0.72993E-01 

0.73852£>01 

0.74605E-01 

•.75242C-01 

0.7S7S7C-01 

0.76152E-01 

0.76431E-01 

0.76694E-01 

0.7S502E-O1 

0.7S060C-01 

0.7SS21E-01 

0.74978E-01 

0.74481E-01 

0.740S5E-01 

0.737t7E-«l 

0.73421E-01 

0.72907E-01 

0.72370C-01 

f.71961E-il 

0.71707C-01 

0.71280E-01 

•.70402£>fl 

0.S9182E-01 

•.680S8C-01 

•.67515E-01 

0.fi7195C-fl 

0.e6425E-01 

•.S5101E-91 

•.6377CC-01 

0.S1503C-01 

•.SStlSE-tl 

•.S4974E-01 

•.S1366E>fS 

•.4N01C-01 

f.4SS90C>0i 

0.41798E>01 

0.36977E-01 

0.32885E-01 

0.29336E-01 

0.25074E-01 

0.17342E-01 

0.13034E-01 

0.99688E-02 

0.78692E-02 

0.66145E-02 

0.58911E-O2 

0.54115E-02 

0.51149E-02 

0.48058E-02 

0.46263E-02 

0.46385E-02 

0.4509OE-02 

O.43464E-02 


0.S2896E-01 

0.52927E-01 

0.S2961E-01 

0.93000E-01 

f.53047C-01 

0.S3099E-ti 

0.53155E-01 

0.53211E-01 

•.53260C-01 

•.53296E<01 

•.S3309£-#l 

0.53292E-01 

0.53236E-01 

0.53134E-01 

0.5298iE-fl 

0.52772E-01 

0.52S07E-01 

0.52188E-01 

0.51812E-01 

0.51389E-01 

0.50924E-01 

0.50423E-O1 

0.49894E>01 

0.48499E-01 

0.47087E-01 

0.45743E-01 

0.44327E-01 

0.43482E-01 

0.42S33E-01 

9.41987E*«1 

0.41525E-01 

0.41208E-01 

0.40810E-fl 

0.40514C-01 

0.40287E-O1 

0.40137E-01 

0.39926E-9I 

0.39482E-01 

0.387S3C-01 

0.38060E-01 

•.37S32E-tl 

t.37181E>tl 

0.386SIE-01 

0.3S064E-01 

0.3S3S1E-01 

0.33917E-01 

0.32284E-01 

0.30569E-01 

0.28741C*#! 

•.S$786C-ft 

0,24O72E'01 

0.21746E-01 

0.20133E-01 

O.18770E-01 

0.1863OE-01 

0.13099E-01 

0.10385E-01 

0.816O9E-02 

0.63752E-02 

0.51095E-02 

O.42095E-02 

0.3569OE-02 

0.30707E-02 

0.27473E-02 

0.26165E-02 

0.24179E-02 

0.22401E*02 

0.217S2E-02 


0.77967E-03 

O.625S1E-03 

0.471S6E-03 

0.31929E-03 

0.16845E-03 

0.2O887E-04 

-0.12492E-03 

•0.26558E-03 

-0.40O79E-03 

-O.52940E-03 

-0.65023E-O3 

-0.7S217E-03 

>0.86421E>83 

-0.95550E-O3 

-0.10354E-02 

>O.11036E>02 

>e.ll599£-02 

-e.l2047E-O2 

-0.12384E>02 

-0.12617E-02 

-0.12756E-02 

-0.12814E-02 

'0.12801E-02 

•0.12548E-02 

-0.12130E-02 

-0.11694E-O2 

-0.11334E-02 

-0.11087E-O2 

-a.l0954E-02 

-0.10908E-02 

-0.10915E-02 

-0.10937E-02 

-0.10941E-O2 

-0.10925E-O2 

-0.nO3lE-O2 

-•.11207E-02 

-0.11203E-02 

-0.10906C-02 

>0.1047OE-02 

-0.10145C-02 

>0.1O033£-02 

-0.99824E-93 

-O.97707E-03 

-0.93808E-03 

“0.897S9E-03 

>f.82378E-03 

-0.71338E-03 

-C.63499E-03 

O.54870E-03 

-0.47859E-03 

>0.41662£-O3 

-O.33758E-03 

-0.2561SE-03 

-0.2OS67E-O3 

-0.16952E-03 

-0.11040E-03 

-0.45149E-04 

-C.17792E-04 

-0.45191E-05 

-0.28847E-06 

0.20O7iE-05 

0.29563E-05 

0.18750E-05 

C.187fi0E~05 

O.18074E-05 

0.15171E-05 

0.19172E-05 

0.18845E-05 

0.14818E-05 


0.69506E-02 

0.69248E'02 

0.89046E-02 

0.68949E-02 

0.68997E-02 

0.69223E-02 

0.69647E-O2 

0.7O275E-02 

0.7U01E-02 

9.72108E-02 

0.732GGE-02 

0.7453GE-O2 

0.75871E-02 

0.77220E-02 

0.78532E-O2 

0.79757E-02 

0.80851E-02 

0.8177SE-02 

0.82504E-02 

0.83022E-02 

0.83324E-O2 

0.83416E-02 

0.S3311E-02 

0.82341E-02 

0.80697E-O2 

0.78775E-02 

0.768G1E-02 

0.75123E-02 

0.73649E-02 

0.72471E-02 

0.71571E-02 

0.70888E-02 

0.69808E-02 

0.68788E-02 

0.68914E-02 

0.67529E-02 

0.6G749E-O2 

0.65137E-02 

0.G2880E-02 

9.6eS35E-02 

0.59G70E-02 

0.58931E-02 

0.57561E-02 

0.55388E-02 

O.53178E-02 

0.4933OE-O8 

O.4408OC-02 

0.39560C-O2 

0.34645C*02 

0.31234E-02 

e.27403E-02 

0.2326SE-02 

O.18402E-02 

0.14S55E-02 

0.12129E-02 

0.90526E-03 

0.47234E-03 

0.27774E-03 

0.16598E-03 

0.10257E-03 

0.69858E~O4 

0.5242SE-04 

0.42022E-04 

0.3S592E-O4 

0.30644E-'04 

0.28248E-04 

0.27362E-04 

0.25350E-O4 

0.23G22E-04 
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RUM N0.» 13 V-STATION- ■ 22 PROBE* 2 VMAX- O.S2®^0E+01 WN WO.- 13 

POSITION MM U UT UT U»U TKE 


0..40343E-01 

0.81983E-ei 

0.1336aC400 

0.ieS26E-f00 

0.20690E-fO0 

0.2485-4E4^00 

0.29018E^00 

«.33182E-t^00 

0.373‘4£E-f00 

0.<415ieE^00 

0.4587-IEtOe 

0.43838E-^00 

0.54O02E>OO 

0.S8166E+00 

8.62330E-fee 

9.66^9<4E+O0 

0.70658E+0O 

0.7-18S3E^00 

0.78986E«00 

0.83iS0E400 

0.873HE400 

0.9H78E40O 

0.95642E+0O 

0.99806E-^00 

9.103S7E-^01 

0.10813E+01 

0.1123OE+01 

0.11646E+01 

0.13063E-»01 

0.12479E401 

0.13895E+O1 

9.13312E401 

0.13728E+91 

0.14145E401 

0.14561E+01 

0.15e03E^01 

0.16S43E+O1 

0.17e84E«01 

0.187S5E401 

0.1976£E«91 

0.30897E4^91 

9.31848E+01 

• .33839Ei^ei 

•.33930£«01 

t.2£912E401 

•.28094E401 

•.30176E4tl 

0.3£25BE4«1 

•.34340E4tl 

t.36422C4fl 


0.779HE+00 
0.77909E->'O0 
0.77905E+00 
0.779O1E4^00 
0.77899E+00 
0.77899E400 
0.77901E^00 
0.77906E4O0 
0.77914E«00 
e.77926E«00 
0.77941E400 
0.77959E400 
0.77982E+00 
0.7S009E+00 
0.78039E400 
O.78074E-»00 
0.78113E400 
0.781S5Ef00 
0.78301E-t-00 
0.78351E>00 
0.78304E400 
0.78361E-f00 
0. 78421 E-i-Oe 
0.78485E+00 
0.78551E+00 
0.78821E+00 
0.78693E-<-00 
0.787£9E-^00 
0.7884SE^00 
0.78929E-t-00 
0.79ei4E^ee 
0.79102E400 
0.79ig2E-fOe 


e.7g28£E-f0e 

8.79383E^e0 

9.79€41E-fO0 

e.79918E-^00 

0.80315E-^O0 

0.80S30E400 

9.80859E-»00 

0.8U97E-^00 

0.81542Em 

0.81890E400 

0.82337E400 

0.8293OE«00 

0.83630E400 

0.843S7E400 

•.8S124E400 

•.85917E400 


•.86706E4ft 

0.87451E400 


0.49674E-01 

0.49670E-01 

0.49£69E'O1 

0.49671E-01 

0.49879E-01 

0.49695E-01 

0.49731E-01 

0.49759E-01 

9.49810E-01 

0.49877E-01 

0.499eiE-01 

0.50062E-01 

9.5ei82E'01 

0.50331E-01 

0.50478E-01 

0.50693E-01 

0.50a45E-01 

0.51052E-01 

0.51273E-01 

0.51505E-01 

0.51745E-01 

0.S1992E-01 

0.52243E-01 

0.52496E-01 

0.527S0E-01 

0.5300SE-O1 

0.53252E-01 

0.53498E-O1 

0.53741E-01 

0.S3981E-01 

0.5421GE-01 

0.54449E-01 


0.54678E-01 

0.54905E-01 

0.55130E-01 

0.55681E-01 

0.56217E>ei 

0.56735E-01 

0.57235E-01 

0.57722E~01 

0.58207E'01 

0.58694E--01 


0.59178E-01 

•.59£35E-01 

0.G0320E-01 

0.G0540E-01 

9.G04S6E-01 

0.60459E-01 

9.6«£33E'01 

•.fi0fi37E-fl 

9.60247C>01 


0.4OS86E401 

0.426S8E401 

•.44750E401 

0.46832E4O1 

0.48914E401 

e. 5O99£E401 
0.5S1G0E«01 
9.59334E>01 
9.G3488E401 
0.G76S2£-^01 
0.71816E^01 
0.759S0E-»01 
0.80144E+01 
0.84308E-»01 
0.88473E>01 

f. 92G3£E401 


0.88138E4O0 
9.88771E400 
0.893G8E400 
9.89943E409 
•.g0492E408 
0.90998E409 
0.91831E+00 
9.g2548£480 
0.93403E«9O 
0.94311E+O0 
0.94931E+00 
0.9S367E+00 
0.9698OE4-00 
0.97023E«Oe 
0.97725E400 
0.98318E4^0e 
a Qoqc:-7F4AA 


0.59701E-61 

e.59282E>01- 

9.5887ee-«l 

O.58209E-01 

0.57313E-01 

0.56309E-01 

0.S3977E-91 

0.53130E-01 

0.49923E-01 

0.46900E'-01 

O.440SGE-01 

0.40716E'-ei 

e.37SI0E-01 

0.33702E-01 

0.29G45E-01 

9.2G732E-01 

a ;>4a<;7r-Ai 


0.36093E-01 

0.3G041E-O1 

0.35991E-01 

0.35944E-01 

0.35900E-01 

0.35SG1E-01 

0.35837E-01 

0.35799E-01 

0.35776E-01 

0.3G759E-01 

0.35746E-01 

0.35738E-01 

9.35734E-01 

0.35732E-01 

0.35731E-01 

O.35732E-01 

0.35733E-01 

0.35734E-01 

0.3573SE-01 

0.35734E-01 

0.35733E-01 

0.35730E-01 

0.3S727E-O1 

0.35723E-01 

0.35719E-01 

0.3571SE-01 

0.35711E-01 

0.35709E-01 

0.35707E-01 

0.3S708E-01 

0.35710E-01 

0.3B71SE-01 

0.35723E-01 

0.35733E-01 

0.35746E-01 

0.35789E-O1 

0.35839E-01 

0.3S882E-01 

0.35903E-01 

0.35894E-01 


0.12123E-03 
0.92523E-04 
0.63987E-94 
0.35782E-04 
0.80415E-05 
-0.19130E-04 
-4.45667E-94 
-a,7154SE-»4 
-0.96774E-04 
'0.12140E-03 
-9.14550E-03 
~9.1G91GE-03 
-0.19246E-03 
-9.21S49E-03 
-0.23832E-03 
-0.26100E-03 
-9.28353E-03 
-0.30592E-03 
-0.3281 IE-03 
-0.35006E-03 
-0.37169E-03 
-0.39291E-03 
-0.413G6E-03 
-0.43386E-03 
-0.45347E-03 
-0.47247E-03 
-0.4908GE-03 
-0.508G8E-03 
-0.52S97E-03 
-0.54279E-03 
-0.5S923E-03 
-3.57535E-03 
-0.59123E-03 
-0.60690E-03 
-0.62240E-03 
-0.G6040E-03 
-0.G9G81E-03 
-0.73067E-03 
-0.76117E-03 
-0.78827E-03 


0.3S8G0E-01 

0.3S820E-01 

0.35797E-01 

0.35S09E-01 

9.35901E-01 

9.358S3E-01 

9.356e0E-01 

0.35227E-01 

0.34872E-91 

0.34748C-01 

0.34387E-01 


-0.81279E-03 

-0.83592E-03 

•0.S5849E-03 

-0.88032E-03 

-9.91599E-03 

-0.93058E-03 

-9.92680E-03 

-0.92238E-03 

-0.92440E-03 

•.92084e-93 

-9.90409E-83 


e.33951E-ei -O.88429E-03 
0.33578E-01 -9.86976E-93 
0.33217E-01 -O.8S3E8E-03 
0.32741E-91 -0.83112E-03 
0.32145E-B1 -0.8OG15E-03 
0.3i504E-»l -0.77907E-03 
0.30277E-01 -0.70304E-03 
0.29200E-01 -0.65093E-03 
0.27993E-01 -0.58367E-03 
0.26387E-01 -0.49852E-03 
0.24747E-01 -0.43G08E-03 
0.233G1E-01 -0.36S37E-03 
0.2i923E-ei -O.30618E-03 
0.20201E-91 -0.24700E-03 
0.18466E-01 -0.19556E-03 
0.17146E-01 -0.15895E-03 
0.15837E-01 -9.12320E-03 


0.37702E-02 

0.37661E-02 

0.37G23E-O2 

0.37591E-02 

0.37568E-02 

0.375SGE-02 

0.375S7E-02 

0.37575E-02 

0.37G10C-02 

0.37GG4E-02 

0.37739E-02 

0.37834E-92 

0.37951E-02 

0.38089C-02 

0.38248E-02 

0.38425E-02 

0.38621E-02 

0.38833E-02 

0.39059E-02 

0.39297E-02 

0.39544E-02 

0.39798E-02 

0.40057E-02 

0.40320E-02 

0.40584E-02 

0.40847E-02 

0.41110E-02 

0.41372E-02 

0.41631E-02 

0.41889E-02 

0.42146E-02 

0.42402E-02 

0.42658E-02 

0.42914E-02 

0.43171E-02 

0.43812E-02 

0.44448E-02 

0.45063E-02 

0.4564SE-02 

0.46202E-02 

0.4G740E-02 

0.47281E-02 

0.47835E-02 

0.48386E-02 

0.49274E-02 

0.49527E-92 

0.49223E-02 

S .489G2E-02 
.48993E-02 
t.48842£-08 
0.48122E-02 
0.47169E-02 
0.46418E-02 
0.45G91E-92 
9.44E03E-92 
0.43181E-92 
0.41G32E-02 
0.38302E-92 
9.35G91E-02 
9.32759E-92 
0.28959E-92 
O.25560E-92 
9.22936E-92 
0.1887GE-92 
0.15439E-02 
0.t2198E-02 
9.10O36E-02 
9.82957E-03 
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